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Description 

The present invention relates to a film thickness controller for use In an extrusion molding apparatus and 
a flowing type molding apparatus such as a film sheet manufacturing apparatus. 
5 A thickness controller is known from US-A-3,557,351 comprising a gauge device and a gauge signal proc- 

essor. The output signal of the gauge signal processor is fed to a controller unit responding to the signal to 
adjust an actuator being coupled to a controlling device incorporated in a material former. The material former 
controls a discharge amount of a material. The gauge is provided for detecting variation thickness of the ma- 
terial after the lapse of a certain time. 
10 From EP-A-0 183 401 the use of a reciprocated gauge is already known. According to FR-A-2,330,058 
state equations are used, while in the publication Chen et al, IEEE vol.3, 20.06.86, pages 1409-1415, the use 
of state equation interferences is explained. 

A conventional film thickness controller is now described briefly. 

An extrusion molding apparatus for manufacturing film or sheet is required to manufacture a molded prod- 

15 uct such as film or sheet having thickness maintained to a predetermined value. An example of a conventional 
apparatus having a die provided with an adjusting mechanism which can adjust thickness of film along the 
width thereof is now described with reference to Figs. 28 to 30. Molten plastic is fed from an extruder 1b (Fig. 
28) to a die 2b. The molten plastic is expanded in a manifold 3b in the width direction perpendicular to paper 
of Fig. 28 showing the die 2b and flows down from a split-shaped outlet 5b of die lips 4b. Then, the molten 

20 plastic flowing dwon from the outlet 5b is cooled by a cooling roller 6b and solidified to be film 7b... so that 
the film is wound on a winder 10b. 

A thickness gauge 11b measures thickness of the film 7b. The thickness gauge 11b utilizes radiation due 
to the natural disintegration of radioactive substance to measure thickness of the film 7b in accordance with 
degree of reduction of the radiation intensity when the radiation passes through the film 7b. The thickness 

25 gauge includes a single detection element which is moved in the reciprocating manner along the width of the 
film 7b to measure thickness of the film 7b along the width. 

It is required that the thickness of the film 7b is maintained to be a predetermined thickness along the width. 
However, since it is difficult that molten plastic passes through a narrow gap of the die 2b in the same speed 
along the width, the thickness of the film 7b is not necessarily identical along the width thereof. 

30 Accordingly, thickness adjusting mechanisms 12b which serve to change a discharge amount of molten 
plastic along the length of the slot of the die lips 4b are disposed dispersedly along the length of the slot of the 
die lips 4b. As the thickness adjusting mechanism 12b, there are the following types, for example: 

(1) Heater type: A multiplicity of heaters are embedded in the die lips 4b along the length of the slot of the 
die lips 4b and are controlled to change a temperature generated therefrom so that the viscosity of the 

35 molten plastic therein is changed and the flowing speed of the molten plastic is changed to control the 
discharge amount of the molten plastic. 

(2) Bolt type: A multiplicity of screws are disposed along the length of the slot of the die lips 4b to change 
a gap space of the discharge outlet 5b of the slot of the die lips 4b mechanically or thermally or electrically 
so that the discharge amount of the molten plastic is controlled. 

40 Accordingly, the thickness of the film 7b can be automatically controlled by adjustment of the thickness 

adjusting mechanism 12b. 

For example, as shown in Figs. 2 and 3, a multiplicity of heaters 12a are embedded in a die 2a at both 

sides of a gap 3a and the heaters 12a are distributed along the width so that the speed of molten plastic flowing 

through the gap 3a is maintained to constant. 
45 At this time, when a temperature of the heater 1 2a which is located in a place where thickness of film 6a 

is thick is reduced, a temperature of molten plastic being in contact with the die 2a is reduced and the viscosity 

of the molten plastic is increased. Accordingly, the flowing speed of the molten plastic therein is reduced. Thus, 

the thickness of a portion of the film 6a corresponding to the place where the temperature of the heater 12a 

is reduced is reduced so that the thicker portion of the film 6a is corrected. Conversely, when the thickness 
50 of the film 6a is small, the temperature of the heater 12 which is disposed in a place where the thickness of 

film is small is increased so that the speed of the molten plastic flowing through the place is increased and 

the thickness of the film 6a therein is increased to correct the thickness of film. 

Fig. 4 is a block diagram of a conventional thickness controller. When a difference between a film thickness 

measured by a thickness gauge 10 and a set value for the film thickness is applied to a controller 1 3, the con- 
55 troller 13 supplies a command to a heater 12a to change a temperature of heat generated by the heater 12a. 

When the temperature of the heater 12a is changed, the flowing speed of the molten plastic in the die 2a is 

changed so that thickness of a portion of the film corresponding to the place where the temperature of the 

heater is changed can be controlled. 
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Fig. 30 is a block diagram of a conventional thickness controller for one operating terminal device of the 
thickness adjusting mechanism 12b. A controller 13b is supplied with a difference between a thickness b of a 
portion of the film measured by the thickness gauge 11b and a set value a of thickness. The controller 13b 
calculates an amount of operation for the adjusting mechanism 12b corresponding to the portion of the film 

5 measured by the thickness gauge 11b and supplies it to the adjusting mechanism 12b. When the mechanism 
1 2b is operated, a discharge amount of molten plastic in the die lips 4b is changed and thickness of the portion 
of the film controlled by the mechanism 12b is changed to effect the thickness control. The thickness control 
over the whole width of the film can be made by provision of the number of the control loop blocks of Fig. 30 
corresponding to the number of places in which the thickness control is performed. 

10 The conventional film thickness controller as described above has drawbacks as follows: 

(1) There is a dead time Li due to movement of the film from the outlet of the die to the thickness gauge 
10 until variation of thickness of the film is detected by the thickness gauge 10 after the variation has been 
produced at the outlet of the die. 

The control system performs a calculation each time a thickness gauge which is reciprocated along the 

15 width of the film reaches an end of the film in which the thickness gauge completes reading of thickness data 
of the film along the width thereof. Consequently, an operation until the thickness gauge reaches the end of 
the film after the thickness gauge has measured thickness of a portion of the film takes a time, which is a 
dead time L 2 until the control system starts the calculation actually after the thickness data has been obtained. 
Accordingly, a dead time after the operation amount for the operating terminal device has been changed and 

20 its influence has been detected as a thickness data until the detected thickness data is employed to perform 
the calculation is a sum of the dead time L, described in (1) and the dead time L 2 described above. 

As described above, the conventional film thickness controller has a drawback of producing a large dead 
time. Description is now made to problems due to these drawbacks. 

25 A. Problem due to large dead time: 

Fig. 5(a) is a block diagram of a thickness control system in the case where the controller of Fig. 4 involves 
the deat times Li and L 2 . Fig. 5(b) is a block diagram of a thickness control system in which the dead times 
are combined to one equivalent time. A general feedback control system does not contain such a dead time, 
30 while the thickness control system contains such a large dead time (Li and L2) as shown in Fig. 5(b). 

Consequently, since there is a large phase delay due to the dead time, a gain of a controller can not be 
increased even if phase compensation is effected in order to attain stability in the control system. Accordingly, 
the high-speed response and the steady-state accuracy of the control system are deteriorated. Further, the 
thickness of the film is always influenced by an external disturbance due to variation of an adjacent die lip 
35 adjusting mechanism. B. Problem due to interference effect: 

In Fig. 30, when an operating terminal device of a portion of the conventional adjusting mechanism 12b 
is operated, the thickness of a portion of the film corresponding to an adjacent operating terminal device is 
changed. Accordingly, the operating terminal device of the portion of the adjusting mechanism and the control 
loop for controlling thickness of a portion of film corresponding to the position of the operating terminal device 
40 interfere with each other. Consequently, the following problems occur 

(1) Even if the stability of the control loop shown in Fig. 30 is ensured, since operation of an operating 
terminal device of the adjusting mechanism 1 2b is influenced by the control loop which controls thickness 
of the film corresponding to an adjacent operating terminal device, the control loops interfere with each 
other and the stability of the whole control system is not ensured when the thickness of the film is controlled 

45 over the whole width of the film. Accordingly, in order to eliminate the influence of the mutual interference, 
the gain of the controller 13b is reduced so that the control system has a low-speed response. 

(2) Conversely, when it is considered to design a stable control system constituting a multi-variable system 
in consideration of the mutual interference between the operating terminal devices of the adjusting mech- 
anism 12b, the control system becomes a very large system since a hundred or more operating terminal 

so devices are usually disposed in the longitudinal direction of the slot of the die lips 4b and there are detected 
values of the film thickness equal to the number of the operating terminal devices. Accordingly, it is difficult 
to design such a large system with ensured stability. 

3. OBJECT AND SUMMARY OF THE INVENTION 

55 

It is an object of the present invention to provide a film thickness controller having a control device which 
solves the problems due to the dead time in a film thickness controller having a large dead time. 
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A. SUMMARY OF FIRST INVENTION 

A film thickness controller for use in an extrusion molding apparatus and a flowing type molding apparatus 
of film including a die having a mechanism which controls a discharge amount of molten plastic along the width 

5 of the film and a thickness gauge for detecting a variation of thickness of the film after the elapse of a dead 
time Li corresponding to a time required for movement of the film between the die and the thickness gauge, 
comprises a subtracter for producing a difference between a thickness value detected by the thickness gauge 
in a predetermined position along the width of the film and a set value of thickness in the predetermined pos- 
ition, an integrator for time-integrating the difference of thickness produced by said subtracter, a memory for 

10 storing past time sequence data of operation amounts of an operating terminal device during a time equal to 
a sum of the dead time and a time L 2 until the thickness gauge reaches an end of the film after detection 
of thickness in the predetermined position, an operational calculator for producing the past time sequence data 
of operation amounts of the operating terminal device stored in said memory and an estimated value of a state 
variable at a time earlier than a time when the set value of the detected thickness value of film has been in- 

15 putted by a dead time L, a state shifter for receiving an output of said integrator and an output of said operational 
calculator and multiplying a coefficient for shifting the state by the dead time L to produce a state estimated 
value at a predetermined time, a state prediction device for receiving the past time sequence data of operation 
amounts of the operating terminal device stored in said memory to produce variation of a state based on es- 
tablishment of input from a certain time to a time after the lapse of the dead time L, an adder for adding an 

20 output of said state shifter and an output of said state prediction device to produce the state estimated value 
at the predetermined time, and an operation amount commander for multiplying a state estimated value at a 
certain time produced from said adder by a state feedback gain to produce an operation amount command 
value of said operating terminal device. 

According to the first invention, a multiplicity of heaters are disposed along the width of the film to control 

25 a temperature of molten plastic which is material of the film, and the thickness gauge detects actual thickness 
of the film at a position downstream of the flowing film and corresponding to the position of the heater in the 
width direction of the film. A difference between the detected actual thickness and a set thickness is calculated 
by the subtracter and is time-integrated by the integrator while a correction command is fed back. Thus, a tem- 
perature of the heater is controlled and a temperature of the molten plastic is controlled to adjust the fluidity 

30 thereof so that thickness of the film is always maintained within the set value. The phase delay due to the dead 
time is corrected by estimation of the past state corresponding to the dead time by the operational calculator 
and time-integration during the time corresponding the past state by the state shifter and the state prediction 
device. 

35 B. SUMMARY OF SECOND INVENTION 

Afilm thickness controller for use in an extrusion molding apparatus and a flowing type molding apparatus 
including a die having a slot along which a plurality of operating terminal devices of a discharge amount ad- 
justing mechanism of molten plastic are disposed and a thickness gauge for detecting variation of thickness 
40 after the lapse of a dead time corresponding to a time required for movement of the film between the die and 
the thickness gauge, comprises a thickness data memory for storing thickness data produced by the thickness 
gauge, a distributor for receiving an output of said thickness data memory and an arrival end identification 
signal which is produced by the thickness gauge to identify whether the thickness gauge reaches either of both 
ends of the film, a plurality of basic control means for receiving an output of said distributor and the arrival end 
45 identification signal produced by the thickness gauge, a plurality of command value memories each receiving 
an output of each of said plurality of basic control means, a superposition adder for receiving an output of each 
of said command value memories, and an operation value memory for receiving an output of said superposition 
adder and for supplying an output of said operation value memory to said basic control means. 
According to the second invention, the following operation is attained. 
50 (1) The thickness gauge measures thickness of the film while moving in the reciprocating manner along 
the width of the film. Since the film is moved at a certain speed, the thickness gauge measures the film 
thickness along a locus as shown in Fig. 27. Accordingly, the thickness gauge produces thickness data of 
a portion of the film corresponding to each operating terminal device sequentially and also produces an 
arrival end identification signal indicating whether the thickness gauge reaches one end (A) or the other 
55 end (B) when the thickness gauge reaches an end of the film. 

(2) The thickness data memory stores thickness data of the film which are measured by the thickness 
gauge over the whole width of the film and which are thickness data of each portion of the film correspond- 
ing to each of the operating terminal devices. 
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(3) The distributor receives the arrival end identification signal of the thickness gauge and further receives 
the thickness data over the whole width of the film from the thickness data memory at the same time as 
receiving of the arrival end identification signal. The distributor supplies a set of predetermined number 
of thickness data from the received thickness data to a predetermined basic control system to be described 

s later. 

(4) Each of basic control systems (control means) receives the set of thickness data supplied from the 
distributor and the arrival end identification signal from the thickness gauge and further receives data set 
from the operation amount memory described later to calculate operation amount command values for a 
plurality of adjacent operating terminal devices containing a predetermined operating terminal device so 

10 that the thickness of a portion of the film corresponding to the predetermined operating terminal device 
is controlled to a predetermined value stably. 

(5) The command value memories store the operation amount command values of the plurality of operating 
terminal devices calculated by the corresponding basic control systems, respectively. 

(6) The superposition adder receives contents of the command value memories storing the operation 
15 amount command values of the basic control systems corresponding to each of operating terminal devices 

and effects superposition, addition and average operation to the command values of each of the operating 
terminal devices to define final command values of each of the operating terminal devices. 

(7) The operation amount memory stores the operation amount command values of each of the operating 
terminal devices defined by the superposition adder retroactively to the past by a time corresponding to 

20 a sum L(=L 1 +L 2 ) of the dead time L, of the thickness gauge and a time L 2 required for movement of the 

thickness gauge from the position corresponding to each of the operating terminal devices to an end of 
the film. 

As described above, the basic control systems can control thickness of the film corresponding to each of 
the heaters (operating terminal devices) containing in the own systems to a predetermined value and can con- 
25 trol thickness over the whole width of the film by combination of the basic control systems. 

4. BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram of a controller according to a first embodiment of the first invention; 
30 Fig. 2 schematically illustrates a configuration of a conventional film manufacturing plant; 
Fig. 3 is a front view showing an conventional arrangement of heaters embedded in a die; 
Fig. 4 is a block diagram of a conventional film thickness controller; 

Fig. 5 is a block diagram of a conventional film thickness controller containing dead time, in which Fig. 
5(a) is a block diagram having separate blocks expressing dead times L, and L 2 , respectively, and Fig. 

35 5(b) is a block diagram having a combined block expressing a sum of the dead times Li and L 2 ; 

Fig. 6 illustrates a correspondence of positions of five heaters and give thickness detection positions; 

Fig. 7 is a block diagram expressing a dynamic mathematical model of film thickness; 

Fig. 8 shows a locus of a thickness gauge for detecting thickness of film; 

Fig. 9 is a diagram illustrating a time interval of calculation and time-integration section; 

40 Figs. 10, 11 and 12 are diagrams illustrating various time- integration sections; 

Figs. 13 and 14 are graphs showing simulation results using an apparatus of the first embodiment of the 
first invention (when a set value of thickness is changed and when external heat is added to a heater, re- 
spectively); 

Figs. 15 to 21 are diagrams concerning a second embodiment of the first invention, in which; 

45 Figs. 15, 16 and 17 are diagram illustrating time intervals of calculation and time-integration sections; 

Figs. 18 and 19 are graphs showing simulation results illustrating effects in the case where an average 
dead time L is used as an integration section of a state shifter and a state prediction device; and 
Figs. 20 and 21 are graphs showing simulation results (when a set value of thicknes is changed and when 
external heat is added to a heater, respectively); 

so Fig. 22 is a block diagram showing a configuration of a controller of a first embodiment of the second in- 
vention; 

Fig. 23 is a block diagram expressing a dynamic mathematical model of a film thickness manufacturing 
process of the first embodiment of the second invention; 

Fig. 24 is a block diagram showing a configuration of a basic control system of the embodiment; 
55 Fig. 25 illustrates an application procedure of the basic control system of Fig. 24 to thickness control points; 

Fig. 26 illustrates a correspondence of positions of five arbitrary operating terminal devices and five thick- 
ness detection positions of the embodiment; 

Fig. 27 illustrates a locus of a thickness gauge which is reciprocated to detect thickness of film in the em- 
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bodiment; 

Fig. 28 schematically illustrates a configuration of a conventional film manufacturing plant; 
Fig. 29 illustrates an arrangement of operating terminal devices embedded in a die of Fig. 28; 
Fig. 30 is a block diagram showing a configuration of the conventional film thickness controller; 
5 Figs. 31(a) to 34(a) are graphs showing simulation results of the embodiment when a set value of thickness 

is changed and Figs. 31(b) to 33(b) are graphs showing simulation results of the embodiment when ex- 
ternal heat is added to a heater; 

Fig. 35 illustrates a discrete time for determining a gain matrix of an operational calculator of the second 
embodiment; and 

10 Figs. 36(a) to 39(a) are graphs showing simulation results of the second embodiment when a set value of 

thickness is changed, and Figs. 36(b) to 39(b) are graphs showing simulation results of the second em- 
bodiment when external heat is added to a heater. 

5. DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

15 

A1. First Embodiment of First Invention 
(a) Transfer Function Matrix G(s) 

20 In order to explain the embodiment, referring to Fig. 6, control of thickness 3' of film to a predetermined 

value is considered by employing five heaters 1 to 5 and thickness values V to 5' of film measured by thickness 
gauges 10 located corresponding to the heaters 1 to 5. The reason that heaters 1, 2 and 3, 4 adjacent to the 
heater 3 are considered in order to control the thickness 3' is to set a control system taking interference of the 
heaters 1, 2 and 4, 5 to the thickness 3' in consideration. Although there are many heaters on both sides of 

25 the heaters 1 and 5, it is considered that influence to the thickness 3' by the heaters disposed outside of the 
heaters 1 and 5 is as small as negligible. Amounts of heat generated by the heaters 1 to 5 are u 1 (t), u 2 (t), u 3 (t), 
u 4 (t) and u 5 (t), respectively, and measured values of thickness V to 5' are y^t), y 2 (t), y 3 (t),y 4 (t) and y 5 (t), re- 
spectively. 

When Laplace transforms of U|(t) and y,(t) (i=1-5) are Uj(s) and Y,(t) (i=1-5), respectively, U,(s) and Y,(s) 
30 are related to each other by the following transfer function matrix G(s): 
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g^s) is a transfer function which introduces temporal variation of thickness 3' when only an amount of heat 
generated by the heater 3, for example, is changed. g 2 (s) is a transfer function which introduces temporal va- 

so riation of thickness 3' when only an amount of heat generated by the heater 2 or 4 is changed. g 3 (s) is a transfer 
function which introduces temporal variation of thickness 3' when only an amount of heat generated by the 
heater 1 or 5 is changed. The equation (1) does not contain a dead time due to movement of the film from an 
outlet of the die to the thickness gauge. Accordingly, g t (s), g 2 (s) and g 3 (s) are expressed by a rational function 
of Laplacian operator s. Further, non-diagonal items of the transfer function matrix G(s) of the equation (1) 

55 express interference to thickness by the adjacent heaters. 
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(b) State Equation 

When the relation between the input U,(s) and the output Y,(s) (1-1-5) of the equation (1) is expressed, the 
following state equation in the canonical form which is convenient for design of the control system is employed: 
5 x(t) = Ax(t) + Bu(t) (2) 

y(t) = Cx(t) (3) 

where x is a state vector, u is an input vector in which u(t)=[ui(t), u 2 (t), u 3 (t), u 4 (t), u 5 (t)F (where T expresses 
transposition), and y is an output vector in which y(t)=[y 1 (t),y 2 (t),y3(t),y 4 (t),y5(t)] T . The state equations (2) and 
(3) are controllable and observable. 
10 The state vector means a vector consisting of a set of variables in which a state of the system is defined 
when the vector is obtained. 

The input vector means a set of variables which is a boundary condition of the state equation for the state 
vector, and the system is controlled by controlling the input vector. 

The output vector means a vector consisting of a set of measured amounts defined by the state vector, 
15 and the system is controlled by measuring the output vector. 

The term "controllable" means that the state vector can be controlled by the input vector. 

The term "observable" means that the state vector can be found by the output vector. 

(c) Output Equation Taking Dead Time into Consideration 

20 

Assuming that a dead time due to movement of the film from an outlet of the die to the thickness gauge 
is Li and a time required for movement of the thickness gauge from the thickness measurement point 3' to an 
end of the film is L 2 , the whole dead time L of the output vector y is given, as shown in Fig. 5(b), by: 

L = Li + L 2 (4) 
25 Accordingly, the output equation (3) is expressed by: 

y(t) = Cx(t-L) (5) 

The relation between the input u(t) (amount of heat generated by the heater) and the output y(t) (detected value 
of the thickness gauge) is shown in Fig. 7 on the basis of the equations (2) and (5). 

30 (d) Arrival End Identification signal 

The thickness gauge measures thickness of film while moving in reciprocating manner along the width of 
the film. Since the film is moved at a certain speed, the thickness gauge measures the film thickness along 
a locus as shown in Fig. 8. When the position of the thickness 3' is shown by ® in Fig. 8, the dead time L 2 
35 due to movement of the thickness gauge is expressed by a time L 2 corresponding to movement between ® 
and® in Fig. 8. 

On the other hand, when the control calculation is made at an end of the film shown by ©of Fig. 8, the 
dead time L 2 due to movement of the thickness gauge is expressed by a time L 2 " corresponding to movement 
between <§T and ©in Fig. 8. As seen from Fig. 8, since the times L 2 ' and L 2 " are different, the control system 
40 which controls the thickness 3' to a predetermined value is characterized in that the dead time L of Fig. 7 in 
the case where the control calculation is made at the nd of the film shown by (B) of Fig. 8 is different from 
that in the case where the control calculation is made at the end of the film shown by (C) of Fig. 8. Accordingly, 
the thickness gauge produces the arrival end identification signal d indicating an end at which the thickness 
gauge arrives. 

45 

(e) Integrator and Output ^(t) thereof 

In order to avoid influence of external disturbance due to thermal conduction from an adjacent heaters to 
control the thickness 3' to a set value, an integrator is introduced to integrate deviation s (t)=r 3 (t)-y 3 (t) between 
50 the detected value y 3 (t) of the thickness 3' and the set value r 3 (t). In the following description, the set value 
r 3 (t)=0. 

The integrator integrates the deviation e (t) until the current time L However, the deviation can be actually 
integrated only until the time (t-L) because of the dead time L. Accordingly, an output X,(t) of the integrator is 
expressed by the following equation: 
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10 



Xi(t) = J*c (t )dr =J^c L (r )dr + J*£ l (t )dr 

Xi(t) = -^C3X(r )dr -J*_C 3 X(r )dr (6) 

where C 3 expresses the third line of C matrix of the equation (3). 
(f) Augmented System State Vector 5f(t) 

15 The first term of the right side of the equation (6) is time-integration of a value capable of being obtained 
from the thickness gauge until time t and accordingly it can be calculated. However, the integrated value of 
the second term of the right side can not be^btained and the time integration can not be calculated as it is. 
Accordingly, in order to obtain prediction of X,(t) at time t, an augmented system as follows in which X,(t) is 
contained in the state variable is considered. 

20 From the equation (6), the following equation is obtained: 



25 



Xi ( t )=-C 3 X( t-L)-C 3 X( t)+C 3 X< t-L)= 
-C 3 X( t) 



From the equations (2) and (7), 
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By using the state vector X(t)=[Xi(t) f X(t)] T of the augmented system, the equation (8) is expressed as follows: 
X( t)=AX( t)+Bu( t) (0) 
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(g) Feedback Gain Matrix 

If the state feedback gain matrix for the equation (9) is F=[fi .FJ, the input u(t) is given by 



55 
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u(t) = -FX(t) = -[fi t F 2 ] 

5 



10 -f ,X, ( t) - F 2 X(t) (11) 



where f t expresses the first column of F matrix. Ifthe feedback gain matrix F is defined so that all characteristic 
values of matrix (A-BF) are in the stable region if Xi(t) and X(t) are obtained, the thickness y 3 (t) can be controlled 
15 to a predetermined value on the basis of the input u(t) stably. Further, since the matricesA and B are not in- 
fluenced by the dead time, this design method can determine the feedback gain matrix F as if it is a system 
having no dead time L and can obtain the high-speed response and the steady-state accuracy of the control 
system. 



Xi ( t) 
X( t) 



20 (h) Calculation of X,(t) and X(t) 



The problem is whether X,(t) and X(t) can be calculated or not If X,(t) and X(t) at the current time t can not 
be obtained, the stable control can not be obtained in the case of the above mentioned feedback gain matrix 
F, and the high-speed response and the steady-state accuracy of the control system are both deteriorated. 
25 The problem (2) in the prior art can not be solved. 

X,(t) and X(t) are obtained as shown in the equation (12) by initializing the time (t-L) and integrating the 
equation (9) from the time (t-L) to the time t Since the input u(t) is already known, the state values X,(t) and 
X(t) are estimated by performing the integration retroactively to the past by the time L. 



30 



Xi ( t) 



Xi(t-L) 



> Bu(t)dt 



35 



X( t) 



X(t~L) 



40 



45 



= e' 



X, ( t-L) 



X( t-L) 



* *t§£* t_T> Bu( x ) d T 



(12 ) 



(i) Calculation of X,(t-L) and X(t-L) 

50 Xj(t-L) of the first term of the right side of the equation (12) is expressed on the basis of the equation (7) 
as follows: 

X, (t-L) = - J*CaX( t )dr (13) 



55 Since the right side of the equation is calculable and is an integrated value of control deviation of the output 
y 3 (t) at the current time t, the equation (1 3) te expressed by: 

XKt-L) = X,(t) (14) 

where X,(t) is an integrated value of control deviation of the detected value y 3 (t) of the thickness 3'. 
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X(t-L) can be estimated as follows: 
From the equations (2) and (5), 

X(t-L) = Ax(t-L) + Bu(t-L) (15) 
y(t) = Cx(t-L) (16) 
5 A variable © (t) defined by the following equation is introduced. 

a>(t) = X(t-L) (17) 
From the equations (15) to (17), the following equations are obtained. 

d>(t) = Aoo (t) + Bu(t-L) (18) 
y(t) = Co(t) (19) 

10 The operational calculator for the equations (18) and (19) is designed to obtain an estimated value X(t-L)=© 
(t) from the detected thickness signal y(t). 

(J) Dead Time L and Calculation Period 

15 Since the calculation is performed each time the thickness gauge reaches the point ® or ©as shown in 

Fig. 8, that is, at regular intervals of time T. 

The relation of the dead time L and the period T of performing the control calculation is described. It is 
assumed that the position ® of the thickness 37 exists near the end ©of the film as shown in Fig. 8. When 
the control calculation is made at the end ® of the film, the whole dead time L of the equation (4) is large 
20 since the dead time L 2 ' is large. On the other hand, when the control calculation is made at the end ©of the 
film, the whole dead time L is small since the dead time L 2 " is small. 

In the embodiment, the dead time L is classified into the following two cases. A case to which the dead 
time L belongs is determined by the arrival end identification signal produced when the thickness gauge reach- 
es the end of the film. 
25 Case1:2T=iL<3T 
Case 2: T < L < 2T 

(k) Discrete Equation 

30 It is necessary to change the equations (18) and (19) to discrete equations for each time interval T and 
design the operational calculator. 

(1)Case1 (2T^L<3T) 

35 As shown in Fig. 9, it is assumed that the control calculation is performed at time V 3 to t*. ^ . It is assumed 
that at the time t^, the output vector y(k+1) is obtained as a set of thickness data and the input vector u(k) is 
maintained constant during time t* to t^. 
From the equation (18), the following equation is derived. 

40 1k x 

<*> Ct K ^v)=e A ( t K -i-t K >w (t K )>+J K *e A (tK*i-T ) 



Bu(r -Ddr <20> 
If the following variable is introduced, the equation (20) is expressed by the equation (21). 



7} = t K ~ 1 - T 



cj (t H ^i) = e AT o> ( t K >+£e^Bu(t K ^i-? -L>d* (21) 

The integration of the right side of the equation (21) means that the double-line portion of Fig. 10 is integrated. 
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Accordingly, the equation (21) is expressed by 

u) (t K -i)=e AT w ( tO+Te^Bd? u(k-2) 



+ S^e^Bdv u(k-3> (22) 



10 The following variables <j> , r 1f r 2 are introduced. 

* - e^ (23) 



r x = ,fe^Bd>7 <24) 



15 



r 2 = .f o e A ?Bd77 C25) 

20 

m = 3t - L (26) 

If the discrete value co (tj is expressed by © (k), the equation (22) is expressed by 

(d (k + 1) = <(> © (k) + r t u(k - 3) + r 2 u(k - 2) (27) 
The discrete equation of the equation (19) is given by 
25 y(k + 1) = C©(k+1) (28) 

By designing the operational calculator for the equations (27) and (28), the estimated value © (k+1) at time 
t=t k+1 is obtained from the following two equations. 



30 u (k+1) = 0 & (k) + T !U(k-3) + r 2 u(k-2) (29) 

u> (k+1) = S (k+1 )+K[y(k+l )-Cw (k+1 )] (30) 

35 where K is a feedback gain matrix of the operational calculator. 

According to the equations (29) and (30), the state © (k+1 9 at time t=t k ^ 1 can be estimated from the set of thick- 
ness data y(k+1) at time t=W The estimated error© (k)=© (k)-© (k) at this time is expressed by 

©(k + 1) = «> - KC<|>)©(k) (31) 
Accordingly, if the gain matrix K of the operational calculator is defined so that all the eigen values of matrix 

40 (<j> -KC<|> ) exist in the stable region, the estimated error can be reduced with the lapse of time. 

(2) Case 2 (T<L<2T) 

As shown in Fig. 10, it is assumed that control calculation is performed at time t k . 2 . and t k+1 . The 
45 integration of the right side of the equation (21) means that the double-line portion of Fig. 10 is integrated. The 
discrete equation at this time is expressed by 

©(k + i) = <|>©(k) + r 1 u(k-2) + r 2 u(k- 1) (31) 
m of the equations (24) and (25) is given by 

m = 2T - L (32) 

so The estimated value © (k+1) at time t=t k ^ 1 is obtained from the following two equations. 



to (k+1) = 0o> (k) + T ,u(k-2) + T 2 u(k-l) (33; 



55 



cj (k+1) = cj (k+1 )+KCy(k+l )-Cw (k+1 ) ] (34) 
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The equation of the estimated error is the same as the equation (31) and the same thing as the case 1 is ap- 
plicable in order to reduce the estimated error with the lapse of time. 

From the foregoing, the estimated value of the state (WL) at t=t k+1 can be obtained in the following se- 
quence. 

5 (1) If time t=tjc +1 is a termination time of the calculation execution period T and it is understood from the 

arrival end identification signal produced from the thickness gauge that the thickness gauge reaches the 
end (B) of the film as shown in Fig. 8, & (k+1) is calculated from the equations (29) and (30) and the es- 
timated value xOm-Dso (k+1) of x(t k+1 -L) is obtained. 

(2) If time t=tjc +1 is a termination time of the control calculation execution period T and it is understood from 
10 the arrival end identification signal produced from the thickness gauge that the thickness gauge reaches 
the end (C) of the film as shown in Fig. 8, © (k+1) is calculated from the equations (33) and (34) and the 
estimated value x(t k+1 -L)=© (k+1) is obtained. 

(1) Calculation of Second Term of Equation (12) 

15 

The final thing to do is to obtain the integration term of the right side of the equation (12), that is, 

1= J*e* c k_r> Bu( t )dr is 

20 

The integration I is to predict variation of the state 



X( t) 



by the input u(t) from time (t-L) to time t. At this time, the dead time L is classified to the following two cases. 
30 A case to which the dead time L belongs is determined by the arrival end identification signal produced from 
the thickness gauge. 

Case1:2T^L<3T 

Case 2:T<L<2T 

35 (1) Case 1.(2T^ LOT) 

In the integration I, the double-line portion of Fig. 11 is integrated. 

xn t w* i ~ L+T-ni- 

I (k+1 ) = $ e A (t K M-T )Bdru(K-2) + 

t K - i-L+T-m+T- 
J e A ( U-i-T )Bd r u(k-l ) + 

45 t k - i -L+T-m 



jt K ~i-L+T-m+2T- ( U ^_ T )BdT u(k) (35) 
tx-i-L+T-m+T 

so 

m = 3T - L 

If the following variable is introduced, the integration I is expressed by the equation (37). 

n = tk ♦ , - t (36) 
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Bd<j u(k-l)+£e**Bda u(k) (37) 

5 o 



(2) Case 2 (T<L<2T) 

In the integration I, the double-line portion of Fig. 12 is integrated. 

10 

U... -L+T-m- 
Kk+l) = ! e A (t K .i-T )Bdr u(k-l) + 

Ik* i-L 



15 

S e A (t K *i-T ) Bd t u(k) 
t K ~i-L+T-m 



20 If the variable r\ of the equation (36) is introduced, the integration I is expressed by 



Kk+l) = e ;T ^"e^Bd(7 u(k-l)+J pBdJ u(k) (38) 

25 



(m) Estimated Value of State Value [X,(t), X(t)] T 



30 From the equations (12), (14), (29), (30), (33), (34), (37) and (38), the estimated value [X x (k+1), X (k+1)] T 
of the state value [X|(t), X(t)P" at current time t=t k+1 is obtained from the following equation. 



X, (k*l) 


= e Xl - 


X, 


(k+1) 


♦ I(k+1) 


(39) 


X(k+1) 




/> 

\ 


(k+1) 







40 

(n) Means for Executing Calculation 

Fig. 1 is a block diagram of a controller implementing the first invention. In the first embodiment, each of 
45 blocks is operated as follows. 

(1) The detected value y(k+1) of film thickness (vector consisting of y^k+1, y 2 (k+1), y 3 (k+1), y 4 (k+1) and 
y5(k+1)) is obtained through a thickness gauge 10 and a sampler 100 at the control calculation execution 
time t=tfc +1 of the time interval T. The sampler 100 closes for each calculation execution time t=i k+u that 
is, the sampler 100 closes each time the thickness gauge 10 reaches the end ® or ©of the film shown 

so in Fig. 8. Further, when the thicknes gauge 1 0 reaches the end ® or © of the film, the gauge 1 0 produces 
the arrival end identification signal d which indicates the end to which the gauge has reached. 

(2) The detected value y 3 (k+1) of the detected film thickness value y(k+1) is supplied to a subtracter 101 
which produces thickness deviation s (k+1)=r 3 (k+1)-y 3 (k+1) between the detected value y 3 (k+1) and a set 
value of thickness r 3 (k+1). 

55 (3) The integrator 1 02 is supplied with the thickness deviation e (k+1 ) from the subtracter 1 01 and produces 
a time-integrated value of the thickness deviation from the following equation. 

X,(k + 1) = X,(k) + 0.5ft ♦ , - t^e (k) + e (k + 1)} (40) 
where e (k) is a thickness deviation at the last thickness detection time (t=t k ) and X<(k) is an output of the 
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integrator 102 at t=tk. 

The integrator 102 includes a function of an external disturbance compensator and serves to com- 
pensate external heat varying the thickness y 3 with heat generated by the heater so that the thickness y 3 
is always maintained to be a set value. 

5 (4) When the thickness gauge reaches either end of the film, the thickness gauge produces the arrival 

end identification signal d. to (k+1) is calculated from the equations (29) and (30) or (33) and (34) in re- 
sponse to the identification signal d. More particularly, in the equation (29) and (30) for the past time se- 
quence data of heat generated by the heater stored a memory 1 04, u(k-3) and u(k-2) are supplied to the 
operational calculator, while in the equations (33) and (34), u(k-2) and u(k-1) together with the detected 

10 film thickness value y(k+1) are supplied to the operational calculator, which produces an estimated value 
X(WL)=£ (k+1 ) of the state variable at time tk +1 -L) earlier than time t^ by the dead time L determined 
by the arrival end identification signal d produced by the thickness gauge. 

(5) In the calculation of the first term of the right side of the equation (39), the state estimated value [X,(k+1 ), 
© (k+1 )] T at time (W L) is multiplied by a coefficient e^ for shifting the state by the dead time L to obtain 

15 the state estimated value e^Xrfk+l), © (k+1)] T at time W That is, the output X,(k+1) of the integrator 
102 and the output © (k+1) of the operational calculator 103 are supplied to state shifter 105, which mul- 
tiplies them by the coefficient for shifting the state by the dead time L determined by the arrival end iden- 
tification signal d of the thickness gauge to obtain the state estimated value at time t k+1 . Since the mag- 
nitude of the dead_time L is different depending on the end of the film which the thickness gauge reaches, 

20 the coefficient e^ is different depending to the position of the thickness gauge upon control calculation 

execution, that is, the arrival end identification signal d of the thickness gauge. 

The state shift by the input u(k) applied in time domain for only the dead time L is expressed by the 
second term l(k+1 )of the right side of the equation (39) and correction therefor is made by a state prediction 
device 106. 

25 (6) The second term !(k+1 ) of the right side of the equation (39) expresses an amount of shift of states for 
time sequence input data u(k-2), u(k-1) and u(k) applied to the time domain from time (WL) to time t* +1 . 
I(k+1) is calculated from the equation (37) or (38) depending on the end of the film which the thickness 
gauge reaches, that is, depending on the arrival end identification signal produced by the thickness gauge. 
More particularly, the past time sequence data of the heat generated by the heater (in this case, three 

30 data of u(k-2), u(k-1) and u(k)) determined by the magnitude of the dead time L stored in the memory 104 
are supplied to the state prediction device 106 and the state variation amount l(k+1) by the input u(k) from 
time (tfrM-LJtotimetk+i. 

(7) Output e^Xrfk+l), © (k+1)] T of the state shifter 105 and output l(k+1) of the state prediction device 

2 A 

35 are added in adder 1 07 which produces the state estimated value [X i (k+1), X(k+1)P" at time t^. Thus, 
although the operational calculator 103 can obtain only the state estimated value at time t k+r L due to the 
dead time L, the state estimated value at time tk +1 can be obtained by integration in the state shifter 105 
and the state prediction device 106 for the dead time L. Influence of phase delay due to the dead time L 
can be eliminated by this operation. 

40 (8) An amount u(k+1 ) of heat generated by the heater from time t k+t to next time tk +2 for control calculation 
is defined by the following equation using state feedback gain (f 1p FJ. 

u(k+1) = -fi^ I (k+1)-F 2 X(k+1) (41) 

A A 

45 The adder 107 supplies the state estimated value [ X (k+1), X(k+1)] T at time t +1 to a commander 108 for 

£ A 

heat generated by the heater. The commander 1 08 multiplies the state estimated value [ x (k+1), X(k+1)] T 
by the state feedback gain to define a command value of heat generated by the heater. 
(9) The above control calculation is executed after the next detected value y(k+2) of film thickness is ob- 
so tained from the sampler 100 at time t=^ 2 of calculation execution when the thickness gauge is moved 
along the width of the film after the time period T and reaches the opposite film end. 



(o) Example of Design 

55 As a first actual example, an example of design is described in the case where transfer functions g^s), 

g 2 (s) and g 3 (s) are given by the following equations: 

, * __ 0.044 ( * 2 \ 

9U ' " S 3 + 2.1S2 + 2.6S + 0.05 1 ' 
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M _ 0.0009 / 43 v 

" s< +2.4S3 + 2.7S2 + 0.25S + 0.0015 1 ; 

n M - 0.00002 m) 

- S5 + 2.4S 4 + 2.8S3 + 0.31S2 + 0.0084S + 0.0004 v ' 
5 u,(t)(i=1-5) is variation (Kcal/s) of heat generated by the heater, and yi(t)(i=1-5) is variation (cm) of thickness 
of film at the position of the thickness gauge corresponding to the position of the heater. The dead time due 
to movement of the film and times L, ' and L 2 " required for movement of the thickness gauge from the thickness 
control point 3' to the film end assume the following values. 
L, = 30 seconds 
10 L 2 ' = 15 seconds 

L 2 " = 1.5 seconds 

It is assumed that the thickness control point 3' exists at the end ©of the film as shown in Fig. 8. The control 
calculation execution period T assumes the following value. 

T = 16.5 seconds (45) 

15 In order to design the control system, it is necessary to express the relation between the input u(t) and the 
output y(t) of the equation (1 ) and obtain the controllable and observable state equations (2) and (3). G(s) con- 
stituted of g^s), g 2 (s) and g 3 (s) of the equations (42) to (44) can be expressed by an equation of the 77th de- 
gree, while the controllable and observable equation has been found to be an equation of the 39th degree. 
Accordingly, the equations (2) and (3) of the 39th degree are obtained from G(s). 



20 



35 



50 



(1) Decision of State Feedback Gain Matrix F 



The state feedback gain matrix F of the equation (11) is obtained as a solution of an optimum regulator 
problem for the state equation (8) extended to the equation of the 40th degree on the basis of the equation 
25 (2). Since the equation (8) is a state equation of a continuous time system, the equation is changed to a discrete 
state function with the sampling period T=16.5 seconds and aregulator solution is applied. A proper estimation 
function is employed to obtain the state feedback gain matrix F and as a result the following values are obtained 
as the eigen values of the matrix (A-BF). 

0.876± 0.02i, 0.79, 0.50± 0.07i, 0.60± 0.09i, 0.60± 0.06i, 0.51 
30 Further, 30 eigen values other than above are not described since the absolute value thereof is less than 0.1 
and attenuation is fast. Since all eigen values are within a circle having a radius of 1, stable control can be 
attained. Since the eigen value having the slowest attenuation is 0.88± 0.02i, the stabilization time Ts can be 
predicted as about 10 minutes from (0876) 35 = 0.01 with the definition of control error 1% as follows. 
Ts = Tx35 = 16.5x35 sec. = 577.5 sec. = 9.6 min. 



(2) Decision of Feedback Gain K of Operational Calculator 



The feedback gain matrix K of the operational calculator of the equation (30) or (34) is obtained for the 
state equation (27) or (31) of the 39th degree and the output equation (28) of the fifth degree. The gain matrix 
40 K is obtained as a solution of the optimum regulator problem so that the matrix T -(C<t> yKJ) has a stable eigen 
value. As a result of obtaining the gain matrix K using a proper estimation function, the following values are 
obtained as eigen values of the matrix ($ -KC<t> ). 

0.9077± 0.0002i, 0.9076, 0.9075, 0.9075, 0.772 ± 0.0001 i, 0.722, 0.722, 0.722, 0.576± 1x10-«i, 0.576± 
1x10-*i, 0.232, 0.232, 0.232, 0.232, 0.232 
45 30 eigen values other than above concentrate to the origin. Since all the values are within a circle having 

a radius of 1 , the estimated error can be reduced with the lapse of time. Since the eigen value having the slow- 
est attenuation is 0.9077, the time To required for attenuation of the estimated error to an initial 1% can be 
predicted from (0.9077) 46 *=? 0.01 as follows. 

To = Tx45 = 16.5x45 sec. = 742.5 sec. = 12.4 min. 



(p) Simulation Example (1) 



Fig. 13 shows an example of simulation result obtained by calculation using the gain matrices F and K 
obtained above. Fig. 1 3(a) shows variations (variations of detected values of the thickness gauge ) of five thick- 
55 ness values y y to y 6 versus time when the set value of thickness y 3 is changed stepwise by 0.02mm. Fig. 1 3(b) 
shows variations of amounts to u 6 of heat generated by the heaters in the same condition as Fig. 13(a). 

Since calculation is made after the execution period of 16.5 seconds of calculation after the set value of 
thickness has been changed, variation of heat generated by the heater occurs after 1 6.5 seconds from change 
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of the set value of thickness. An amount of heat generated by the heater is maintained to the same value until 
1 6.5 seconds elapse and the next calculation is made. The calculation is made on the basis of a newly detected 
value of thickness after 16.5 seconds to change an amount of heat generated by the heater. Accordingly, an 
amount of heat generated by the heater changes stepwise as shown in Fig. 1 3(b). 

5 On the other hand, variation of the detected thickness value is detected after the lapse of the dead time 

L of 31.5 seconds after the amount of heat generated by the heater has been changed after the lapse of 16.5 
seconds from the change of the set value. That is, variation of thickness is detected after the lapse of 16.5 + 
31.5 = 48 seconds after the set value of thickness has been changed. Thickness y 3 is exactly changed to a 
set value and the change is substantially symmetrical to the thickness y 3 . Variation of heat generated by the 

10 heater u 3 is largest, variations by the heaters u 2 and U4 are largest next to the heater u 3 , and variations of the 
heaters and u 5 are smallest 

The stabilization time is about 1 8.5 minutes which is considerably large as compared with the stabilization 
time of 12.4 minutes calculated by the eigen value of the operational calculator (the stabilization time by the 
eigen value of the regulator is stilt shorter). This is based on the following reason. 

15 In order to prevent the command value of heat generated by the heater from being changed largely for 
each calculation, the command value is defined with weight added as follows. 

u d ,k = Wu d , k . , + (1-W)u k (46) 

where 

u dfk : command value of heat defined by the calculation time t=t kf 

20 u dfk . y command value of heat defined by the last calculation time tFt*. 1f 

u k : command value of heat calculated at the calculation time t=tk, and 

W: weight coefficient. 

In this simulation, W=0.8. This means that when the calculation period T=16.5 seconds is considered, a 
time delay corresponding to a delay of first order having a time constant of 74.65 seconds is added to the heat 

25 commander. Accordingly, it is considered that the stabilization time of thickness control of Fig. 13 is larger than 
the stabilization time estimated by the eigen value of the operational calculator. Then, even if the thickness 
control is in the stabilization state, the command value of heat changes for each calculation. The reason is 
because the magnitude of the dead time L of the first term e^ of the right side of the equation (39) which is 
one of the calculation equations is different in one end ® and the other end © of the film as shown in Fig. 8. 

30 When the present control system is applied actually, the same calculation equation as that applied to the 
thickness y 3 is applied to each of thicknesses y 1p y 2 , y 4 , and y 6 and each command value of heat may be pro- 
duced as a sum of results of the calculation equations. It will be understood that the control system considerably 
reduces influence of the dead time since the time required for stabilizing variation of thickness when heat gen- 
erated by the heater is changed stepwise without control of heat is about 10 minutes. 

35 

(q) Simulation Example (2) 

The second actual example is now described with reference to Fig. 14, which shows a control result when 
external heat of 8.4 wattage is applied to the heater U3. Fig. 14(a) shows variations of thickness values y 1 to 

40 y 5 versus time, and Fig. 1 4(b) shows variations of heat Uj to u 5 generated by the heaters versus time. As shown 
in Fig. 14(a), although the thickness y 3 is once increased by the external heat of the heater u 3 , the thickness 
y 3 is returned to the original set value by changing the amounts of heat generated by the heaters u, to u 5 and 
the stabilization time is about 18.5 minutes in the same manner as Fig. 13. It is understood that variation due 
to the external disturbance is exactly compensated by introducing the integrator in the present control system. 

45 The thickness values y 2 and y 4 are once increased by influence of external heat through thermal conduction 
along the width of the die. The thickness values y t and y 5 are also influenced similarly, although the influence 
is small as compared with y 2 and y 4 . In order to cancel the influence of such external heat, reduction of heat 
generated by the heater u 3 is largest, reduction by the heaters u 2 and u 4 is largest next to the heater u 3 , and 
reduction by the heaters and u 6 is smallest. When external heat is applied to the heater u 3 , other thickness 

50 values y 1f y 2 , y 4 and y 5 are also changed, although such interference effect can be canceled by applying the 
same calculation equation as that for the thickness value y 3 to each of the thickness values y 1f y 2 , y 4 and y 5 . 

A2. Second Embodiment of First Invention 

55 (a) Relation to First Embodiment of First Invention 

In the second embodiment, the process that the same elements as in the first embodiment are utilized, 
the equations (1 ) to (1 9) are derived, the operational calculator for the equations (1 8) and (1 9) is designed and 
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the estimated value X(t-L)=© ((t) of X(t-L) is obtained from the detected value of thickness y(t) is quite identical 
with that of the first embodiment. 

In the second embodiment, a known © (y obtained by the calculation performed in the step just before 
the current step is used to obtain © (t^). 

5 

(b) Dead Time 

In the second embodiment, the calculation is also executed each time the thickness gauge reaches the 
end ® or ©of the film as shown in Fig. 8. That is, the calculation is executed at regular intervals of time T. 
10 The time T is a time required for movement of the thickness gauge along the width of the film. 

On the other hand, the dead time L for the position ® of thickness 3 is different depending on the end ® 
and the end ©. That is, 

Dead time by calculation for the end ® : L B =L,+L 2 " (47) 
Dead time by calculation for the end © : L C =L 1 +L 2 " 
is It is apparent that Lb>I_c for the position ® . In the description below, it is assumed that L B and L c satisfy 
T<L<2T. 

(cj Known © (tk) 

20 Fig. 1 5 is a diagram for explaining the calculation for the end ® for obtaining the estimated value © (W 

L B ). It is assumed that the calculation is made at time V2 to tfc* 1a Further, it is assumed that the input u(k) of 
the heater is maintained to constant from tk to W 

It is assumed that the thickness gauge reaches the end ® at time t k+1 . Accordingly, it is considered that 
the thickness gauge has reached the end ©at the past time ^ before time t* +1 by time T. Thus, it is assumed 

25 that the estimated value © (t k )= X(t k -Lc) has been obtained in the calculation for the end ©executed at time 
t k . In the calculation for the end ® executed at time t^ u the known © (tj is employed to obtain the estimated 
value © (t kH )=X(t k+1 -L B ). A 

Fig. 16 is a diagram for explaining the calculation for the end ©for obtaining the estimated value X(t k+1 - 
l_c). It is assumed that the thickness gauge reaches the end ©at time tk +1 and the thickness gauge has reached 

30 the end ® at the past time tk earlier than time t^ by time T. In the calculation for the end © executed at time 
t k+1 , the known © (W is employed to obtain the estimated value &(t k+1 )=X(t k+1 -L c ). 

As seen from Figs. 15 and 16, the time interval T of the calculation is constant, while since the dead times 
L B and Lc for the ends @ and ©, respectively, are different, the time difference expressing the estimated val- 
ues ©(tj and co(tk +1 ) is different from the time interval T. Accordingly, the estimated value © (t^) is obtained 

35 from the equations (18) and (19) as follows. 

(d) Discrete Equation (1 8) for the End ® 

A A 

The calculation for the end ® shown in Fig. 1 5 is first considered. The known estimated value ©(tk)=X(tk- 
40 L c ) is expressed by X(to). The estimated value © (VO^WLb) to be obtained is expressed by X(t0. The es- 
timated value X(t,) of state variable at time tj is estimated from the equation (18) on the basis of the known 
X(to) and inputs after time to. 



X(t x ) = e A < tl - lo> X(to)+J^e* ctl - T5 Bu(T )dr (48) 



to = ^ - t t = t^ + , - L B , t! - to = T - L B + Lc 
When a new variable i\ = t r x is introduced, the estimated value x(*i) is transformed as follows: 



55 XUx) = e A< * l - %0, X<to>*J*e^§u(ti-ir )dt (49) 

Since the integration of the right side of the equation (49) means that the double-line portion of Fig. 15 is in- 
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tegrated, the equation (49) is expressed by the following equation. 



s X(t x ) = e A (T-L B *Lc)X(to)+ J 0 e A, Bd7 u(k-l) 

T-Lb+Lc 

*S e^BdT? u(k-2) (50) 

JPb 

10 

m B = 2T - L B 

where u(k-1) is a heater input from time t^ t to t*. and u(k-2) is a heater input from time %. 2 to V v If x"<(ti)= 
£f (tk +1 ) and X(to)-© (W, the equation(50) is expressed by 

15 

<J (k+1) = e A (T-L B -L c >w (k)+ J™% A "Bd* u(k-l) 

20 

T *~ L ♦ L. 

+ J e A71 Bd?? u(k-2) (51) 

IDb 

25 where the discrete value © (y is expressed by © (k). 

(e) Discrete Equations (18) for the End © 

The calculation for the end ©shown in Fig. 16 is considered. At this time, in the equation (49) 
30 to = tk - Lb, t t = t k + ! - l_c, t t - to = T - Lc + L B 

By integrating the double-line portion of Fig. 16, the following equation is obtained. 

X(t x ) = e A (T-L c +L B )X( t 0 >+J o e A7l Bd7 u(k-l> 

35 

+ S T ~ L t'^Bdv u<k-2> (52) 

JJ)c 

40 m c = 2T - l_c 

If x (*i)= u> and X(to)=A (tk). the estimated value ©(ta) is given from the equation (52) by 

45 u (k+1) = e A (T-L B -L c )u> (k)+I o e A7f Bd? u(k-l) 

T-Lb+Lg 

+ J e A7I Bd?7 u(k-2) (53) 

so 

(f) Discrete Equation for Equation (19) 



55 



The discrete equation for the equation (19) is given by 

y(k+1) = C©(k + 1) (54) 
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(g) Calculation of Estimated value © (k+1) 

By designing the operational calculator in accordance with the equations (51 ), (53) and (54), the estimated 
value © (t k+1 ) at t=t k+1 is obtained as follows. 
5 The calculation equation of the estimated value © (k+1) for the end ® is given by 

u> (k+1) = 0 b<*> (k) + S 0 e^hdy u(k-l) 



10 

+ , 

IDs 



+ J T ^ L ?^Bd? u(k-2) (55) 



15 u (k+1) = c5 (k+1) + K B Cy(k+l)-C< A ) (k+1)] (56) 

where $ q-q^T-Lq^-Lc) 

K B =gain matrix of the operational calculator. 
The calculation equation of the estimated value © (k+1 ) for the end © is given by 



20 



25 



30 



u (k+1) = 0 c w (k) + J Q e A7? Bd? u(k-l) 



T-Lo+Lb 

+ S e* n hdy u(k-2) (57) 

IDc 

(j (k+l) = Z (k+l) + K c Cy(k+l ) -Cg> (k+1)] (58) 



where <t> c-e^T-Lc+U) 

Kc=gain matrix of the operational calculator. 
According to the equations (55) and (58), the state variable © (k+1) at t=tk +1 can be estimated by a set of 
35 thickness data y(k+1) at t=tfc +1 . 

(h) Estimated Error © (k) 

At this time, the estimated error © (k)=© (k)- © (k) is expressed by the following equation: 
The estimated error © (k+1) for the end ® is given by 

©(k+1) = [<J>b~ KbC* b]o (k) (59) 
The estimated error © (k+1) for j_he end © is given by 

©(k+1) = [(t>c-KcC<t>c]a)(k) (60) 
Accordingly, if the gain matrices Kb and Kc of the operational calculator are defined so that all eigen values 
of the matrices b -KbC<|> b ] and c-KcC<j> c ] are in the stable domain, the estimated value can be reduced 
with the lapse of time. 

(i) Summary of Calculation of Estimated Value of © (k+1) 

so From the foregoing, the estimated value of the state X^-L) at t=tK +1 can be obtained in accordance with 
the following sequence. 

(1) When t=tK +1 is a termination time of the period T of calculation execution and it is discriminated on the 
basis of the arrival end identification signal produced by the thickness gauge that the thickness gauge 
has reached th end ® ofthe film shown in Fig. 8, © (k+1) is calculated from the equations (55) and (56) 

55 and the estimated value X(tk +1 -Le)=© (k+1) of XOWU) is obtained. 

(2) When t=tn +1 is a termination time of the period T of calculation execution and it is discriminated on the 
basis of the arrival end identification signal produced by the thickness gauge that the thickness gauge 
has reached th end ©ofthe film shown in Fig. 8, © (k+1) is calculated from the equations (57) and (58) 
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and the estimated value Xfo+t-Lc) 55 © (k+1) is obtained. Thus, the first term of the right side of the equation 
(12) can be calculated. 

(j) Integration of Second Term of Equation (12) 

The final thing to do is to obtain the integration term of the right side of the equation (12), that is, 



10 



■ ■ Ce 



AC t - X 5 



Bu( t )dr 



This integration term I is to predict variation of the state 



15 



Xi (t) 
X( t) 



20 by the input u(t) from time (t-L) to time t. The integration I is to integrate the double-line portion of Fig. 17. 



t K * x -L+T-m 
I(k+1> = 5 e A (t K ~i-r )Bdr u(k-l) 



25 



♦ J tHt e A (t K .i-T ) Bd t u(k> 
t*^i-L+T-m 



30 m = 2T - L 

If a new variable tj =tk +1 -T is introduced, the integration l(k+1) is expressed by 



I (k+1 ) = e* T I L e" Af Bdau(k-l) 



35 



* J T e***Bda u(k) 

o 



(61) 



Since the dead time L is different depending on the calculation for the ends ® and ©, the equation (61) is 
40 described as follows. 

The integration l B (k+1) for the end ® is given by 



45 



l»(k*l) = e AT W tf Bdau(k-l 



+ £e x **Bda u(k) 



(62) 



so The integration l c (k+1 ) for the end © is given by 



I c (k+1 ) = e AT jI; c 3 T Bdau(k-l) 



55 



+ ^e^BdJ u(k) 

o 



(62) 
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(k) Estimated value [X i (k+1), X (k+1)] T 

From the equations (12), (14), (55), (56), (57), (58), (62) and (63), the estimated value [Xi (k+1), X (k+1)] T 
of the state value (Xrft), X(t)Y at the current time t=t* +1 is obtained by the following equations. 
The estimated value for the end ® is given by 



10 



15 



Xi (k+1) 



X(k+1) 



Xi(k+1) 



u> (k+1) 

The estimated value for the end © is given by 



+ I B (k+l> 



(64) 



20 



25 



Xi (k+1 ) 



X(k+1) 



Xi (k+1) 



cj (k+1) 



+ I c (k+1) 



(64) 



(1) Discontinuity of Estimated Value [X x (k+1 ), X (k+1 )]T 



30 



35 



40 



45 



50 



When the calculation equation for the estimated value [X r (k+1), X(k+1 )] T of the state value at time t^ is 
changed depending on the end ® or ©as described in the equation (64) and (65), since the dead times L B 
and Lc are different and change stepwise, the estimated value is not continuous each time the equation is 
changed. When the dead time La is larger than the dead time L c , the estimated value of the equation (64) is 
larger than that of the equation (65) and accordingly the estimated value by the equation (11) repeatedly 
changes in the pulse manner. Fig. 18 shows a simulation result when the estimated value is calculated using 
the equations (64) and (65). In Fig. 18, the time interval of calculation t=22.5 seconds, the dead time L B =39 
seconds and L c =37 seconds. Fig. 1 8 shows a control result when external heat of 8W is applied to the heater 
u 3 stepwise. Fig. 18(a) shows variations of thickness values y, to y 6 versus time, and Fig. 1 8(b) shows variations 
of amounts u 1 to u 5 of heat generated by the heaters versus time. As shown in Fig. 18(a), the thickness y 3 is 
once increased by the external heat, although the thickness y 3 is returned to the original set value by changing 
the amounts of heat generated by the heaters to u 6 . However, heat generated by the heaters is repeatedly 
changed in the steady state and the thickness is also slightly changed repeatedly. When the position of thick- 
ness y 3 approaches the end of the film, since a difference between the dead times I_b and Lc is increased, a 
width of variation of heat generated by the heater is increased and variation of thickness is also larger when 
the estimated equations (64) and (65) are employed. 

(m) Average Value L of Dead Time 

In order to improve the above drawback, an average value of Lg and Lc, that is, L=(Le+Lc)/2 is employed 
as the dead time used in the equations (64) and (65). The equation of the estimated can be used in common 
for the ends ©and © 
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Xi(k+1) 



X(k*l) 



Xi(k+1) 



u> (k+1) 



+ I(k*l> 



(66) 



10 



_- L-T - r _ 
l(k+l) = e AT J e A ' T Ddff u(k-l) 



f5 



+ re A ' r Bd(r »(k) 

O 



(67) 



L = (L B + Lc)/2 (68) 



20 (n) Simulation Example 



25 



Fig. 19 shows a simulation result when the equations (66), (67) an (68) are used as the equation of the 
estimated value with the same condition as in Fig. 18. Variation in the steady state of heat generated by the 
heater is eliminated. 

(o) Means for Executing Calculation 



Fig. 1 is a block diagram of a controller implementing the first invention. In the second embodiment, each 
of blocks is operated as follows. 

30 (1) The detected value y(k+1) of film thickness (vector consisting of y^k+1, y 2 (k+1), y 3 (k+1), y 4 (k+1) and 
y 5 (k+1)) is obtained through the thickness gauge 10 and the sampler 100 at the calculation execution time 
t=^ +1 of the time interval T. The sampler 100 closes for each calculation execution time t=t k + 1f that is, the 
sampler 100 closes each time the thickness gauge 10 reaches the end ® or ©of the film shown in Fig. 
8. Further, when the thickness gauge 10 reaches the end ® or ©of the film, the gauge 10 produces the 

35 arrival end identification signal d which indicates the end which the gauge has reached. 

(2) The detected value y 3 (k+1) of the detected film thickness value y(k+1) is supplied to a subtracter 101 
which produces thickness deviation s (k+1)=r 3 (k+1)-y 3 (k+1) between the detected value y 3 (k+1) and a set 
value of thickness r 3 (k+1). 

(3) The integrator 102 is supplied with the thickness deviation e (k+1 ) from the subtracter 1 01 and produces 
40 a time-integrated value of the thickness deviation from the following equation. 

X,(k + 1) = X,(k) + 0.5ft ♦ i - W{e (k) + e (k + 1)} (69) 
where s (k) is thickness deviation at the last thickness detection time (t=tjj and X,(k) is an output of the 
integrator 102 at t=t k . 

The integrator 102 includes a function of an external disturbance compensator and serves to com- 
45 pensate external heat varying the thickness y 3 with heat generated by the heater so that the thickness y 3 
is always maintained to be a set value. 

(4) When the thickness gauge reaches either end of the film, the thickness gauge produces the arrival 
end identification signal d. & (k+1) is calculated from the equations (55) and (56) or (57) and (58) in re- 
sponse to the identification signal d. More particularly, the past time sequence data u(k-2) and u(k-1) of 

so heat generated by the heater stored a memory 1 04 together with the detected film thickness value y(k+1 ) 
are supplied to the operational calculator, which produces an estimated value X (WL)= & (k+1) of the 
state variable at time t( k+1 -L) earlier than time t^ by the dead time L determined by the arrival end iden- 
tification signal d produced by the thickness gauge. 

(5) In the calculation of the first term of the right side of the equation (66), the state estimated value [X|(k+1 ), 
55 © (k+1)] T at time (t k+1 -L) is multiplied by a coefficient e*»- for shifting the state by the average dead time 

L defined by the equation (68) to obtain the state estimated value e^Xrfk+l), © (k+1)] T at time t^. That 
is, the output X,(k+1 ) of the integrator 1 02 and the output © (k+1 ) of the operational calculator 1 03 are sup- 
plied to state shifter 1 05, which multiplies them by the coefficient for shifting the state by the average dead 
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time L to obtain the state estimated value at time The magnitude of the dead time L adopts the average 
value of the dead times for both ends of the film as described by the equation (68). 

The state shift by the input u(k) applied in time domain for only the average dead time L is expressed 
by the second term l(k+1) of the right side of the equation (66) and correction therefor is made by a state 
5 prediction device 106. 

(6) The second term f(k+1 ) of the right side of the equation (66) expresses an amount of shift of states for 
time sequence input data u(k-1) and u(k) applied to the time domain of the average dead time from time 
(t k+1 -L) to time t k+1 . F(k+1) is calculated from the equation (67) using the average dead time L. More par- 
ticularly, the past time sequence data of the heat generated_by the heater (in this case, two data of u(k- 

10 1) and u(k)) determined by the magnitude of the dead time L stored in the memory 104 are supplied to 
the state prediction device 106 and the state variation amount I (k+1) by the input u(k) from time (t k+1 -L) to 
time t k+1 . 

(7) Output e AL [X,(k+1), © (k+1)F of the state shifter 105 and output T(k+1) of the state prediction device 

15 are added in adder 107 which produces the state estimated value [Xi (k+1), X(k+1)P" at time t^. Thus, 

although the operational calculator 103 can obtain only the state estimated value at time t k+1 -L due to the 
dead time L, the state estimated value at time t* +1 can be obtained by integration in the state shifter 105 
and the state prediction device 106 for the dead time L Influence of phase delay due to the dead time L 
can be eliminated by this operation. 

20 (8) An amount u(k+1) of heat generated by the heater from time t k+1 to next time t k+2 for calculation is defined 
by the following equation using state feedback gain (f tl FJ. 

u(k+1) = - ^(k+IJ-F.Xfk+l) (41) 

A A 

25 The adder 107 supplies the state estimated value [ x (k+1), X(k+1)F at time t^ to a commander 108 for 

A A 

heat generated by the heater. The commander 1 08 multiplies the state estimated value [ x (k+1), X(k+1 )] T 
by the state feedback gain to define a command value of heat generated by the heater. 
(9) The above control calculation is executed after the next detected value y(k+2) of film thickness is ob- 
w tained from the sampler 100 at time t=^ 2 of calculation execution when the thickness gauge is moved 
along the width of the film after the time period T and reaches the opposite film end. 



(p) Example of Design 

35 As a first actual example, an example of design is described in the case where transfer functions g^s), 
g 2 (s) and g 3 (s) are given by the following equations: 

0.14 

9l(S) = S3 + 5.5S*+12.5S + 0.25 (69) 



40 



45 



55 



a / s x _ 0.003 (?Q) 

921 ' " S* + 6.4S3 + 1 3.2S2 + 1 .3S + 0.009 v 9 

( v _ 0.00005 , 71 * 

93V S J " s 5 + 6.3^ + 1 3.8S3 + 1.6S2 + 0.04S + 0.0002 1 ' 



u,(t)(i=1-5) is variation (watt) of heat generated by the heater, and y,(t)(i=1-5) is variation (micron) of thickness 
of film at the position of the thickness gauge corresponding to the position of the heater. The dead time U due 
to movement of the film and times IV and L 2 " required for movement of the thickness gauge from the thickness 
control point 3' to the film end assume the following values. 
U = 26 seconds 
L 2 ' = 17 seconds 
L 2 " = 7.5 seconds 
^ Accordingly 

L B = 43 seconds 
Lc = 33.5 seconds 

It is assumed that the thickness control point 3' exists at the end ©of the film as shown in Fig. 8. The control 
calculation execution period T assumes the following value. 

T = 22.5 seconds (72) 

In order to design the control system, it is necessary to express the relation between the input u(t) and the 
output y(t) of the equation (1 ) and obtain the controllable and observable state equations (2) and (3). G(s) con- 
stituted of g 1 (s), g 2 (8) and g 3 (s) of the equations (69) to (71) can be expressed by an equation of the 77th de- 
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gree, while the controllable and observable equation has been found to be an equation of the 29th degree. 
Accordingly, the equations (2) and (3) of the 29th degree are obtained from G(s). 

(1) Decision of State Feedback Gain Matrix F 
5 _ 

The state feedback gain matrix F of the equation (11) is obtained as a solution of an optimum regulator 
problem for the state equation (8) extended to the equation of the 30th degree on the basis of the equation 

(2) . Since the equation (8) is a state equation of a continuous time system, the equation is changed to a discrete 
state function with the sampling period T=22.5 seconds and aregulator solution is applied. A proper estimation 

10 function is employed to obtain the state feedback gain matrix F and as a result the followingvakies are obtained 
as main values for determining a response of control as the eigen values of the matrix (A-BF). 

0.856, 0.8119,0.7755, 0.7618 
Further, other eigen values except above are not described since the absolute value thereof is small and at- 
tenuation is fast Since all eigen values are within a circle having a radius of 1, stable control can be attained. 
15 Since the eigen value having the slowest attenuation is 0.8560, the stabilization time Ts can be predicted as 
about 11 minutes from (0876) 30 = 0.01 with definition of control error 1% as follows. 

Ts = Tx35 = 22x30 sec. = 675 sec. = 11.3 min. 

(2) Decision of Feedback Gain K of Operational Calculator 

20 

The feedback gain matrix K of the operational calculator of the equation (56) or (58) is obtained for the state 
equation (55) or (58) of the 29th degree and the output equation (54) of the fifth degree. The gain matrix K is obtaind 
as a solution of the optimum regulator problem so that the matrices B T -(C<|> bV^b 1 ) and fcfc 7 -^ cY^c 1 ) have a 
stable eigen value. For example, the discrete time of the state equation (55) defining the gain matrix K B is 
25 (T-Lb+Lc)=(22.5-43+33.5)=13 seconds. As a result of obtaining the gain matrix K using a proper estimation 
function, the following values are obtained as main values for determining convergence of the operational cal- 
culator as eigen values of the matrix (<)> B -KaC<|> B ). 
0.9183, 0.9183, 0.9183, 0.9183, 0.9183, 
0.7654, 0.7654, 0.7654, 0.7654, 0.7654, 
30 Other eigen values except above are not described since the absolute values are small and convergence 
is fast. Since all the values are within a circle having a radius of 1, the estimated error can be reduced with 
the lapse of time. Since the eigen value having the slowest attenuation is 0.9183, the time To required for at- 
tenuation of the estimated error to an initial 1% can be predicted from (0.91 SZ) 55 ^ 0.01 as follows. 
T 0 = (T-L B +L c )x55 = 13x55 sec. = 715 sec. = 12 min. 
35 The gain matrix Kc having the stabilization time To of 12 minutes is obtained for the matrix <t> c . 

(q) Simulation Example 1 

Fig. 20 shows an example of simulation result obtained by calculation using the gain matrices F, K B and 

40 Kc obtained above. Fig. 20(a) shows variations (variations of detected values of the thickness gauge ) of five 
thickness values to y 5 versus time when the set value of thickness y 3 is changed stepwise by 5 micron. Fig. 
20(b) shows variations of amounts to u 5 of heat generated by the heaters in the same condition as Fig. 1 3(a). 

Since calculation is made after the execution period of 22.5 seconds of calculation after the set value of 
thickness has been changed, variation of heat generated by the heater occurs after 22.5 seconds from change 

45 of the set value of thickness. An amount of heat generated by the heater is maintained to the same value until 
22.5 seconds elapse and the next calculation is made. The calculation is made on the basis of a newly detected 
value of thickness after 22.5 seconds to change an amount of heat generated by the heater. Accordingly, an 
amount of heat generated by the heater changes stepwise as shown in Fig. 20(b). 

On the other hand, variation of the detected thickness value is detected after the lapse of the dead time 

50 L of 33.5 seconds after the amount of heat generated by the heater has been changed after the lapse of 22.5 
seconds from the change of the set value. That is, variation of thickness is detected after the lapse of 22.5 + 
33.5 = 56 seconds after the set value of thickness has been changed. Thickness y 3 is exactly changed to a 
set value and the change is substantially symmetrical to the thickness y 3 . Variation of heat generated by the 
heater u 3 is largest, variations by the heaters ^ and u 5 are largest next to the heater u 3 , and variations of the 

55 heaters u 2 and u 4 are smallest. This reason is because interference of the heaters u 2 and u 4 to thickness_y 3 is 
reduced. The stabilization time which is estimated by the eigen value determined by the gain matrices F, Kb 
and Kc and is 12 minutes is supported by Fig. 20. There is no variation in heat generated by the heater at the 
steady state, since the equation (66) is employed to compensate the dead time instead of the equation (64) 
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and (65). 

When the present control system is applied actually, the same calculation equation as that applied to the 
thickness y 3 is applied to each of thicknesses y tl y 2 , y 4 , and y 6 and each command value of heat may be pro- 
duced as a sum of results of the calculation equations. It will be understood thatthe control system considerably 
5 reduces influence of the dead time since the time required for stabilizing variation of thickness when heat gen- 
erated by the heater is changed stepwise without control of heat is about 13 minutes. 

(r) Simulation Example 2 

10 The second actual example is now described with reference to Fig. 21 , which shows a control result when 
external heat of 8 watts is applied to the heater u 3 . Fig. 21(a) shows variations of thickness values y 1 to y 6 
versus time, and Fig. 21(b) shows variations of heat u, to u 6 generated by the heaters versus time. As shown 
in Fig. 21(a), although the thickness y 3 is once increased by the external heat of the heater u 3 , the thickness 
y 3 is returned to the original set value by changing the amounts of heat generated by the heaters to u 5 and 

15 the stabilization time is about 1 2 minutes in the same manner as Fig. 20. It is understood that variation due to 
the external disturbance is exactly compensated by introducing the integrator in the present control system. 

The thickness values y 2 and y 4 are once increased by influence of external heat through thermal conduction 
along the width of the die. The thickness values y, and y 5 are also influenced similarly, although the influence 
is small as compared with y 2 and y 4 . In order to cancel the influence of such external heat, reduction of heat 

20 generated by the heater u 3 is largest, reduction by the heaters Uj and u 6 is largest next to the heater u 3( and 
reduction by the heaters u 2 and u 4 is smallest. This is because the reduction in the heaters u 2 and u 4 does not 
influence thickness y 3 so much. When external heat is applied to the heater u 3 , other thickness values y t , y 2 , 
y 4 and y 5 are also changed, although such interference effect can be canceled by applying the same calculation 
equation as that for the thickness value y 3 to each of the thickness values y 1f y 2 , y 4 and y 5 . 

25 

A3. Effects of the Invention 

The present invention is configured as described above and accordingly has the following effects. The in- 
tegrator which time-integrates a difference between a detected value of thickness of film at a predetermined 

30 position and a set value of thickness is introduced and an output of the integrator is fed back to compensate 
an amount of heat generated by the heater for external heat influencing thickness of the film so that thickness 
of the film can be always identical with the set value. Further, in order to avoid large phase delay due to the 
dead time, the state estimated value at time t-L earlier than the current time t by the dead time L is obtained 
by the operational calculator and the state estimated value at time t-L is time-integrated by the state shifter 

35 and the state prediction device during the dead time L so that the state estimated value at the current time t 
can be obtained to remove deterioration of control performance due to the dead time. 



Claims 

40 

1. A film thickness controller which employs state equations to control an extrusion molding apparatus or 
a flowing type molding apparatus of film including a die having a mechanism which controls a discharge 
amount of molten plastic along the width of the film and a thickness gauge for detecting variation of thick- 
ness of the film after the lapse of a dead time L 1? corresponding to a time required for movement of the 

45 film between the die and the thickness gauge, comprising a subtracter (101) for producing a difference 

between a thickness value detected by the thickness gauge (10) in a predetermined position along the 
width of the film and a set value of thickness in the predetermined position, whereat said thickness gauge 
is reciprocated along the width of the film, characterized in that an integrator (102) for time-integrating 
the difference of thickness produced by said substractor (1), a memory (104) for storing past time se- 

50 quence data of operation amounts of a thickness adjusting device (12b) during a dead time L equal to 

the sum of the dead time L% and a time L 2 until the thickness gauge (1 0) reaches an end of the film after 
detection of thickness in the predetermined position, an operational calculator (1 03) for producing the past 
time sequence data of operation amounts of the thickness adjusting device (12b) stored in said memory 
and an estimated value of state variable at a time earlier than the time when the set value of the detected 

55 thickness value of film has been inputted by the dead time L, a state shifter (105) for receiving an output 

of said integrator (1 02) and an output of said operational calculator (1 03) and multiplying each of said out- 
puts of said integrator (102) and operational calculator (103) by a coefficient for shifitng the state by the 
dead time L to produce a state estimated value at a predetermined time, a state prediction device (106) 
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for receiving the past time sequence data of operation amounts of the thickness adjusting device (12b) 
stored in said memory (104) to produce state variation based on establishment of input from a certain 
time to a time after the lapse of the dead time L, and adder (1 07) for adding an output of said state shifter 
(105) and an output of said state prediction device (106) to produce the state estimated value at the pre- 
determined time, and an operation amount commander (108) for multiplying a state estimated value at a 
certain time produced from said adder by a state feedback gain to produce an operation amount command 
value of said thickness adjusting device (12b). 

A controller according to Claim 1, wherein said state equations are given by 

x(t) = Ax(t) + Bu(t) (2) 
y(t) = Cx(t-L) (4) 

where x is a state vector, u is an input vector in which uftHu^t), u 2 (t), u 3 (t), u 4 (t) ? u 5 (t)] T (where T represents 
transposition), y is an output vector in which y(t)=[yi(t). y 2 (t), y 3 (t), y 4 (t), y 5 (t)] T , and the state equations 
(2) and (3) are controllable and observable. 

A controller according to Claim 2, wherein said integrator (102) performs the following calculation: 



7, ct> = - J*cVx<t )dr - ^CJXCT )dr (6) 
where Ci is the i-th row of C matrix. 

A controller according to Claim 3, wherein when <o (t)=x(t-L), said operational calculator (103) calculates 
an estimated value X(t-L) for x(t-L) on the basis of the following equations for q (t): 

cb(t) = A© (t) + Bu(t-L) (18) 
y(t) = Co(t) (19) 

A controller according to Claim 4, wherein said equation (18) is the following discrete equation: 



co (t K ~i> = e A (t k M-t K )w <t K ) I e A ( t K ~i- T ) 

tic 



Bu( T - L) d T 

A controller according to Claim 6, wherein when Ti=t k+1 -t, an estimated value & (k+1) for a> (t^) is calcu- 
lated by the following equations: in the case of 2^L< T, 

tj (k+J) = 0 Z <k) + r xU(k-3)*r 3 u(k-2> (29) 
w (k+l) = £) (k*l )+KCy(k*l)-CtJ Ck+1)] (30) 

where K is a feedback gain matrix of said operational calculator (108), and 
in the case of T<L< 2T, 

cj (k+l) = 0 co Ck) + T x u(k-2)*r a u(k-l) 
£ Ck+l) = u (k+1 ) + KCy(k+l)-Cc3 (k+1)] 

where 

4 = e AT (23) 

r x = J% A ^Bd? (24) 
T 2 = S^e^Bdv (25) 



(33) 
(34) 
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m = 3T - L (26). 

A controller according to Claim 1 , wherein said controller comprises the following constituents: 

(1) a detected value y(k+1) of film thickness (vector consisting of y^k+1), y 2 (k+1), y 3 (k+1), y 4 (k+1) and 
ys(k+1)) is obtained through the thickness gauge (10) and a sampler (100) at calculation execution time 
t=tfc. i of the time interval T; the sampler (100) closes for each calculation execution time t=tk_ 1( that is, 
the sampler (100) closes each time the thickness gauge 10 reaches the end ® or ©of the film shown 
in Fig. 8; further, when the thickness gauge (10) reaches the end @ or ©of the film, the gauge (10) 
produces the arrival end identification signal d which indicates the end which the thickness gauge (10) 
has reached; 

(2) a value y 3 (k+1) of the detected film thickness value y(k+1) is supplied to subtracter (1 01) which pro- 
duces thickness deviation e (k+1)=r 3 (k+1)-y 3 (k+1) between the detected value y 3 (k+1) and a set value 
of thickness r 3 (k+1); 

(3) the integrator (102) is supplied with the thickness deviation e (k+1) from the subtracter (101) and 
produces a time-integrated value of the thickness deviation from the following equation; 

X,(k+1) = XrfkJ + O^. i-t k ){e(k)+6(k+1)} (40) 
where 8 (k) is thickness deviation at the last thickness detection time (t=y and X,(k) is an output of the 
integrator (102) at t=tn; 

the integrator (102) includes a function of an external disturbance compensator and serves to 
compensate external heat varying the thickness y 3 with heat generated by the heater so that the thick- 
ness y 3 is always maintained to be a set value; 

(4) when the thickness gauge reaches either end of the film, the thickness gauge produces the arrival 
end identification signal d; © (k+1) is calculated from the equations (29) and (30) or (33) and (34) in 
response to the identification signal d; more particularly, in the equations (29) and (30) for the past 
time sequence data of heat generated by the heater (108) stored in memory (104), u(k-3) and u(k-2) 
are supplied to the operational calculator (103), while in the equations (33) and (34, u(k-2) and u(k-1) 
together with the detected film thickness value y(k+1) are supplied to the operational calculator (103), 
which produces an estimated value X(t k _ r L)=© (k+1 ) of the state variable at time t( k . r L) earlier than 
time tti by the dead time L determined by the arrival end identification signal d produced by the thick- 
ness gauge; 

(5) in the calculation of the first term of the right side of the equation (39), the state estimated value 
[X,(k+1), © (k+1)] T at time (tk_ r L) is multiplied by a coefficient e*- for shifting the state by the dead 
time L to obtain the state estimated value e^XKk+l), © (k+1)] T at time t* that is, the output X,(k+1) 
of the integrator (1 02) and the output © (k+1 ) of the operational calculator (1 03) are supplied to state 
shifter (105), which multiplies them by the coefficient for shifting the state by the dead time L deter- 
mined by the arrival end identification signal d of the thickness gauge to obtain the state estimated 
value at time t k . ,; since the magnitude of the dead time_L is different depending on the end of the film 
which the thickness gauge reaches, the coefficient e^ is different depending on the position of the 
thickness gauge upon calculation execution, that is, the arrival end identification signal d of the thick- 
ness gauge ; 

the state shift by the input u(k) applied in time domain for only the dead time L is expressed by 
the second term l(k+1) of the right side of the equation (39) and correction therefor is made by state 
prediction device (106); 

(6) the second term l(k+1) of the right side of the equation (39) expresses an amount of shift of states 
for time sequence input data u(k-2), u(k-1) and u(k) applied to the time domain from time (t k . r L) to 
time ^ t ; I (k+1) is calculated from the equation (37) or (38) depending on the end of the film which 
the thickness gauge reaches, that is, depending on the arrival end identification signal produced by 
the thickness gauge, more particularly, the past time-sequential data of the heat generated by the hea- 
ter (in this case, three data of u(k-2), u(k-1) and u(k)) determined by the magnitude of the dead time L 
stored in the memory (104) are supplied to the state prediction device (106) and the state variation 
amount l(k+1) by the input u(k) from time ft. i-L) to time t*. t is produced; 

(7) output 6*4X^+1), © (k+1)] T of the state shifter (1 05) and output l(k+1) of the state prediction device 

2 A 

(106) are added in adder (107) which produces the state estimated value [Xi (k+1), X(k+1)] T at time 
tK.,; thus, although the operational calculator (103) can obtain only the state estimated value at time 
t k _ r L due to the dead time L, the state estimated value at time t k ^ can be obtained by integration in 
the state shifter (1 05) and the state prediction device (1 06) for the dead time L; influence of phase delay 
due to the dead time L can be eliminated by this operation; 
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(8) an amount u(k+1) of heat generated by the heater (12a) from time t*. ^ to next calculation time t^ 
is defined by the following equation using state feedback gain (f 1t F 2 ) ; 

u(k+1) = -fa^k+IJ-F^k + l) (41) 

the adder (1 07) supplies the state estimated value [ x (k+1 ), X(k+1 )] T at time 1 to a commander (1 08) 

A 

for heat generated by the heater; the commander 108 multiplies the state estimated value [ x (k+1), 
X(k+1)^ by the state feedback gain to determine a command value of heat generated by the heater; 
10 and 

(9) the above control calculation is executed after the next detected value y(k+2) of film thickness is 
obtained from the sampler (100) at time t=tic. 2 of calculation execution when the thickness gauge is 
moved along the width of the film after the time period T and reaches the opposite film end; where at 
equations (29), (30), (33), (34), (37), (38), (39) are given as follows: 

15 £ (k+l)=tf w (k) + T & u(k-3) + r a utk-2> (29) 



u (k+l)=u (k+1 )+K[y(k*l )-Cu (k+l)J (30) 

20 



u<kM> = 0u<k) + r & u<k-2) + r«uCk-L) (33) 

25 

<*> (k+l)=u (k+1 )+K[y(k+l )-Cw (k*l)l (34) 



KkM) = c*<= T ' i^e^fida u(k-2)+e* T J g**" 

Md<r ii (k - 1 ) + f^e^Bda u(k) (37) 

35 

Mk+t) = e^ T J^I^Bda u(k-l) + J g^Bda u(k) (38) 



40 





X t (k+1) 






+ I < k* 1 > 


(39) 


45 


X(k+1 ) 




u> (k + 1 ) 







8. A controller according to Claim 7, wherein the feedback matrix F of said operation value commander (108) 
so is selected so that all eigen values of matrix (A-BF) are in a stable region. 

9. A controller according to Claim 4, wherein said operational calculator (1 08) calculates estimated value © 
(k+1 ) for © (tfc_ 0 from the following equations: 

the calculation equation of the estimated value © (k+1) for the end ® is given by 

55 

cj (k+1) = 0 b cj (k) + i m e A,, Bd7? u(k-l) 

o 
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T-Lb+Lc 

* S e^Bdv u(k-2) (55) 



u (k+1) = u (k-M) + K B Cy(k-»-l)-Cu (k+l)3 (56) 

10 where 

4> B = eA<T-LB + U:) 
K B =gain matrix of the operational calculator (103), 

the calculation equation of the estimated value © (k+1) for the end © is given by 

w (k+1) = <f> c(*> (k) + J o e M Bd? u(k-l) 

+ J^I^Bd? u(k-2> (B7) 



20 



25 



£ (k+1) = u (k*l) * KcCy(k+l)-C<7 (k+l)3 (58) 



where 

*c = e A Cr-Lc + L B ) 
Kc=gain matrix of the operational calculator. 

30 1 0. A controller according to Claim 1 , wherein said controller comprises the following constituents: 

(1) the detected value y(k+1) of film thickness (vector consisting of y^k+1, y 2 (k+1), y 3 (k+1), y 4 (k+1) 
and y 6 (k+1)) is obtained through the thickness gauge (10) and sampler (100) at the calculation execu- 
tion time t=V ! of the time interval T; the sampler (100) doses for each calculation execution time \=\*. 1f 
that is, the sampler (100) closes each time the thickness gauge (10) reaches the en ® or © of the 

35 film shown in Fig. 8; when the thickness gauge (10) reaches the end ® or © of the film, the gauge 

(10) produces the arrival end identification signal d which indicates the end which the gauge has 
reached; 

(2) a value y 3 (k+1) of the detected film thickness value y(k+1) is supplied to subtracter 101 which pro- 
duces thickness deviation e (k+1)=r 3 (k+1)-y 3 (k+1) between the detected value y 3 (k+1) and a set value 

40 of thickness r 3 (k+1); 

(3) the integrator (102) is supplied with the thickness deviation e (k+1) from the subtracter (101) and 
produces a time-integrated value of the thickness deviation from the following equation; 

X,(k + 1) = X,(k) + 0.5(t k . , - 10 {e (k) + s (k + 1)} (69) 
where e (k) is thickness deviation at the last thickness detection time (t=tk) and X,(k) is an output of the 
45 integrator (1 02) at t=t*; 

the integrator (102) includes a function of an external disturbance compensator and serves to 
compensate external heat varying the thickness y 3 with heat generated by the heater so that the thick- 
ness y 3 is always maintained to be a set value; 

(4) when the thickness gauge reaches either end of the film, the thickness gauge produces the arrival 
50 end identification signal d. & (k+1) is calculated from the equations (55) and (56) or (57) and (58) in 

response to the identification signal d; the past time sequence data u(k-2) and u(k-1) of heat generated 
by the heater stored a memory (1 04) together with the detected film thickness value y(k+1) are supplied 
to the operational calculator, which produces an estimated value X(t k . r L)=© (k+1) of the state variable 
at time t( k . r L) earlier than time k , by the dead time L determined by the arrival end identification signal 
55 d produced by the thickness gauge; 

(5) in the calculation of the first term of the right side of the equation (66), the state estimated value 
[X,(k+1), & (k+1)] T at time (t* r L) is multiplied by a coefficient e^" for shifting the state by the average 
dead time C defined by the equation (68) to obtain the state estimated value 6^X^+1), & (k+1)] T at 
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time t k _ the output X^k+1) of the integrator (102) and the output © (k+1) of the operational calculator 
(103) are supplied to state shifter (105), which multiplies them by the coefficient for shifting the state 
by the average dead time E to obtain the state estimated value at time k the magnitude of the dead 
time E adopts the average value of the dead times for both ends of the film as described by the equation 
(68); 

the state shift by the input u(k) applied in time domain for only the average dead time E is ex- 
pressed by the second term T(k+1 ) of the right side of the equation (66) and correction therefor is made 
by state prediction device 106; 

(6) the second term T(k+1) of the right side of the equation (66) expresses an amount of shift of states 
for time sequence input data u(k-1) and u(k) applied to the time domain of the average dead time E 
from time (t k _ r E) to time t ; T(k+1) is calculated from the equation (67) using the average dead time 
E; the past time sequence data of the heat generated by the heater (in this case, two data of u(k-1) 
and u(k)) determined by the magnitude of the dead time E stored in the memory (104) are supplied to 
the state prediction device (106) and the state variation amount l(k+1) by the input u(k) from time (t r E) 
to time t i, 

(7) output e^X^k+l), © (k+1)] T of the state shifter 105 and output l(k+1) of the state prediction device 

are added in adder 107 which produces the state estimated value [X r (k+1), X(k+1)] T at time t k . t ; al- 
though the operational calculator (1 03) can obtain only the state estimated value at time t k . r L due to 
the dead time L; the state estimated value at time t*, , can be obtained by integration in the state shifter 
(105) and the state prediction device (106) for the dead time E; influence of phase delay due to the 
dead time L can be eliminated by this operation; 

(8) an amount u(k+1) of heat generated by the heater from time t^ , to next calculation time t^ is de- 
fined by the following equation using state feedback gain (f 1p F 2 ); 

u(k+1) = -fiXi(k+1)-F 2 X(k + 1) (41) 

a A 

the adder (107) supplies the state estimated value [ X (k+1), X(k+1)] T at time tfc. t to a commander 108 

A 

for heat generated by the heater; the commander (108) multiplies the state estimated value [ x (k+1), 
X(k+1)P by the state feedback gain to define a command value of heat generated by the heater; and 

(9) the above calculation is executed after the next detected value y(k+2) of film thickness is obtained 
from the sampler (1 00) at time t=tn. 2 of calculation execution when the thickness gauge is moved along 
the width of the film after the time period T and reaches the opposite film end, where at equations (55), 
(56), (57), (58), (66), (67), (68) are given as follows : 



Ck*l> = <t> *<3 (k) ♦ J^t^Bd* u(k-l) + J T ~ L ?^Bd? u(k-2> (55) 

uck+n = £<k+i> ♦ KB(y<k*i>-cJ<k*i>] (50) 
where 

^eACT-Le + Lc) 
Ka=gain matrix of the operational calculator 

u (KM) = 0 (k) • /"V'^U u<k-!> 

♦ J T ' L S^Bd» u(k-2> 

w (kM) = Z (kH) * KcCy<k+l )-Cw (k*l)3 
where 

*c = o a (T-Lc + Lb) 
30 



(fi7) 
(58) 
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Kc=gain matrix of the operational calculator. 




X.(k*l) 



e 



♦ I(k*l) 



(66) 



X(k*l) 




♦ 5 T e^ r Bd a u(k) 



(G7) 



L = (Lb + Lc)/2 (68) 



Patentanspruche 
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Vorrichtung zur Kontrolle der Dicke eines Filmes, welche Zustandsgleichungen zur Steuerung eines 
Extrusions-SpritzgieBgerates Oder eines PreBspritzgerStesfur einen Film verwendet, mit einer GuBform, 
welche einen Mechanismus aufweist, der die Ausf luBmenge an geschmolzenem Plastik entlang der Film- 
breite steuert, und einer Dicke-MeBeinrichtung zur Erfassung der Filmdickenabweichung nach Verstrei- 
chen einer Totzeit die einer fur die Bewegung des Filmes zwischen der GuBform und der Dicke- 
MeBeinrichtung erforderlichen Zeit enspricht, und mit einer Subtrahierschaltung (101) zur Erzeugung ei- 
ner Differenz zwischen einem durch die Dicke-MeBeinrichtung (10) in einer vorbestimmten Position ent- 
lang der Filmbreite erfaBten Dicke-Wert und einem Dicke-Soll-Wert in der vorbestimmten Position, wobei 
die Dicke-MeBeinrichtung entlang der Filmbreite hin- und herbewegt wird, 
dadurch gekennzeichnet, 

daB sie einen Integrator (102) zur Zeitintegration der durch die Subtrahierschaltung (101) erzeugten 
Dicken-Differenz aufweist, 

einen Speicher (104) zur Speicherung vergangener Datensequenzen der Betriebsmengen einer Dicken- 
Einstelleinrichtung (1 2b) wahrend einer Totzeit L, welche gleich der Summe ist aus der Totzeit U und einer 
Zeit L 2 , bis die Dicke-MeBeinrichtung (10) ein Ende des Filmes nach dem Erfassen der Dicke in der vor- 
bestimmten Position erreicht, 

einen Betriebsrechner (103) zur Erzeugung der vergangenen Datensequenz der Betriebsmengen der 
Dicke-Einstelleinrichtung (12b), welche in dem Speicher gespeichert sind, 

und eines Schatzwertes der Zustandsvariablen zu einem Zeitpunkt, weicher vor dem Zeitpunkt liegt, an 
dem der Sollwert des erfassten Filmdickewertes mit der Totzeit L eingegeben wurde, 
einen Zustandswechsler (105) fur den Empfang eines Ausgangsignals des Integrators (102) und eines 
Ausgangssignals des Betriebsrechners (103) zur Multiplikation der beiden Ausgangssignale des Integra- 
tors (1 02) und des Betriebsrechners (1 03) mit einem Koeff izienten zum Zustandswechsel bei der Totzeit 
L zur Erzeugung eines Zustands-Schatzwertes bei einer vorbestimmten Zeit, 

eine Zustandsvorhersage-Einrichtung (106) fur den Empfang der vergangenen Datensequenz der Be- 
triebsmengen der Dikkeneinstell-Einrichtung (12b), welche in dem Speicher (104) gespeichert ist, zur Er- 
zeugung einer Zustandsabweichung basierend auf der Ermittlung des Eingangs ab einem gewissen Zeit- 
punkt bis hin zu einem Zeitpunkt nach Verstreichen der Totzeit L, 

einer Addiereinrichtung (107) zur Addition des AusgangsdesZustandswechslers (105) und des Ausgangs 
der Zustandsvorhersage-Einrichtung (106) zur Erzeugung des Zustands-Schatzwertes zur vorbestimm- 
ten Zeit, 

und eine Betriebsmengen-Befehlseinrichtung (108) zur Multiplikation des von der Addiereinrichtung zu 
einem bestimmten Zeitpunkt erzeugten Zustands-Schatzwertes mit einer Zustands-Ruckkopplungs-Ver- 
stSrkung zur Erzeugung eines Betriebsmengen-Befehlwertes fur die Dickeneinstell-Einrichtung (12b). 
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Vorrichtung nach Anspruch 1 bei der die Zustandsgleichungen gegeben sind durch 

x(t) = Ax(t) + Bu(t) (2) 
y(t) = Cx(t-L) (4) 

wobei x ein Zustandsvektor, u ein Eingangsvektor ist mit u(t) = [u^t), u 2 (t), u 3 (t), u 4 (t), u 5 (t)] T (wobei T fur 
Transposition stent), y ein Ausgansvektor ist mit y(t) = [y^t), y 2 (t), y 3 (t), y 4 (t) f y 5 (t)] T , und die Zustandsglei- 
chung (2) und (3) steuerbar und beobachtbar sind. 

Vorrichtung nach Anspruch 3 bei der der Integrator (102) folgende Berechnungen ausfuhrt: 



TiCt) » - J*CiX<T )dT - J^CJXCT )dT (6) 
wobei Cj die i-Reihe der C Matrix ist. 

A 

Vorrichtung nach Anspruch 3 bei der wenn <o(t)= x(t-L), der Betriebsrechner (103) einen SchStzwert X(t- 
L) fur x(t-L) auf der Basis der folgenden Gleichungen fur co(t) berechnet: 

<b(t) ■ Aco(t) + Bu(t-L) (18) 
y(t) = Cco(t) (19) 

Vorrichtung nach Anspruch 4 bei der die Gleichung (18) die folgende diskrete Gleichung ist: 
W (t K -i) = e A ( t K -x-tK>u> <t K > ♦ S e A ( tn-i-T ) 



Bu( r -Ddr 

Vorrichtung nach Anspruch 6 bei der, ri=tk +1 -T wenn ein SchStzwert © (k+1) f ur © (t^) durch die folgenden 
Gleichungen berechnet wird: 
Im Fall da& 2^L<T, 

Z <k*] ) = 0 u (k)>T x uCk-3)*r aUCk-2> (29) 



u <k*l) = <j (k*l)*KCy(k*l)-Cu (k*l>] ( 30 ) 



wobei K eine Ruckkopplungsverstarkungs- Matrix des Betriebsrechners (108) ist, und fur den Fall, dafi 
T<L<2T 

Z <k*l)=0 u (k)>T x u(k-2)*r a u(k-l> (33) 
u <k*l)=u <k*l>*Kty(k*l)-Co (k+l)l (34) 

wobei 

4 = eAT (23) 
r x = jV^Bd? < 24) 



= £e"«Bd* (25) . 
m = 3T - L (26) 

Vorrichtung nach Anspruch 1, bei der die Kontrollvorrichtung folgende Bestandteiie aufweist: 
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(1) ein erfafiter Wert Y(k+1) der Filmdicke (Vektor der bestehend aus y^k+1), y 2 (k+1), y 3 (k+1), y 4 (k+1) 
und y 5 (k+1)) wird erhalten durch die Dicken-MeBeinrichtung (10) und einen Ab taster (100) bei der Be- 
rechnungsausf uhrungszeitt^ des Zeitintervalls T; der Abtaster (100) schlieBt fur jede Berechnungs- 
ausfuhrungszeitt=t k+1 , d.h. der Abtaster (100) schlieBt jedes Mai, wenn die DickemeBeinrichtung (10) 

5 das Ende ©Oder ©des Filmes, wie in Figur 8 gezeigt, erreicht; ferner, wenn die DickemeBeinrichtung 

(1 0) das Ende ® oder ©des Filmes erreicht, erzeugt die DickemeBeinrichtung (1 0) das Ende - erreicht 
-Idendifikationssignal d, welches das Ende, das die DickemeBeinrichtung (10) erreicht hat, anzeigt: 

(2) ein Wert y 3 (k+1) des erfaBten Filmdickewertes y(k+1) wird einer Subtrahierschaltung (101) zuge- 
fuhrt, welchedie Dickenabweichung e(k+1) = r 3 (k+1) - y 3 (k+1) zwischen dem erfaBten Wert y 3 (k+1) und 

10 einem Sollwert der Dicke r 3 (k+1 ) liefert; 

(3) der Integrator (102) erhalt die Dickenabweichung e (k+1) von der Subtrahierschaltung (101) und 
erzeugt einen zeitintegrierten Wert der Dickenabweichung aus der folgenden Gleichung: 

X ( (k + 1) = X,(k) + 0.5(tK . 1 - W {e (k) + e (k + 1)} (40) 
wobei e (k) die Dickenabweichung zum letzten Dicke nerfaBungszeitpunkt (t=tk) und X^k) der Ausgang 
15 des Integrators (102) bei t=^ ist; 

Der Integrator (102) umfaBt eine Funktion fur den Ausgleich externer StQrungen und dient zum Aus- 
gleich externer Hitze, welche die Dicke y 3 verSndert, mit Hitze, welche durch den Erhitzer erzeugt wird, 
so daB die Dicke y 3 immer bei einem Sollwert gehalten wird; 

(4) wenn die DickemeBeinrichtung ein Ende des Filmes erreicht, erzeugt die DickemeBeinrichtung das 
20 Ende-erreicht-ldendif ikationssignal d; ©(k+1) wird aus den Gieichungen (29) und (30) oder (33) und 

(34) in Reaktion auf das Idendifikationssignal d berechnet; genauer gesagt wird in den Gieichungen 
(29) und (30) fur die vergangene Datensequenz der Hitze, welche durch den Erhitzer erzeugt und in 
dem Speicher (1 04) gespeichert sind, u(k-3) und u(k-2) dem Betriebsrechner (103) zugef uhrt, wahrend 
in den Gieichungen (33) und (34) u(k-2) und u(k-1) zusammen mit dem erfaBten Filmdickewert y(k+1) 
25 dem Betriebsrechner (103) zugef uhrt werden, welcher einen SchStzwert X(t k+1 - L)=©(k+1) der Zu- 

standsvariable (tn +1 - L) erzeugt, welcher vor dem Zeitpunkt t^ urn die Totzeit L liegt, welche durch 
das Ende-erreicht-ldendif ikationssignal d bestimmt wird, daB durch die DickemeBeinrichtung erzeugt 
wird; 

(5) bei der Berechnung des ersten Terms der rechten Seite der GleichungJ39) wird der Zustandsschatz- 
30 wert [Xt(k+1), co(k+1)] T zum Zeitpunkt (t k+1 -L) mit einem Koeffizienten e^ zum Wechsel des Zustandes 

bei derTotzeit L multipliziert, urn einen Zustandsschatzwert zum Zeitpunkt t* +1 zu erhalten; d.h. der Aus- 
gang X|(k+1) des Integrators (102) und der Ausgange ©(k+1) des Betriebsrechners (103) werden dem 
Zustandswechsler (105) zugef uhrt, welcher sie mit dem Koeffizienten multipliziert zum Wechsel des 
Zustands nach der Todzeit L, welche durch des Ende-erreicht-ldendifikationssignal d der 
35 DickemeBeinrichtung bestimmt wird, urn den Zustandsschatzwert zum Zeitpunkt t* +1 zu erhalten; da 

die GroBe derTotzeit L in Abhangigkeit des von der DickemeBeinrichtung erreichten Filmendes unter- 
schiedlich ist, ist der Koeff izient e*- in Abhangigkeit der Position der DickemeBeinrichtung bei der Be- 
rechnungsausfuhrung unterschiedlich, d.h. in Abhingigkeit von dem Ende-erreicht-ldendifikationssi- 
gnal d der MeBeinrichtung; 

40 der Zustandswechsel durch den Eingang u(k), welcher in der Zeitdomine nur fur die Totzeit L angelegt 

wird, wird durch den zweiten Term I (k+1 ) der rechten Seite der Gleichung (39) ausgedruckt und die Kor- 
rektur wird durch die Zustandsvorhersageeinrichtung (106) durchgefuhrt; 

(6) der zweite Term l(k+1 ) der rechten Seite der Gleichung (39) druckt einen Betrag des Zustandswech- 
selsfurdieZeiteingabedatensequenzu(k-2),u(k-1)und u(k) aus, die fur das Zeitintervall vom Zeitpunkt 

45 (WL) bis zum Zeitpunkt t^ anliegen; l(k+1) wird aus der Gleichung (37) oder (38) in Abhangigkeit 

des Filmendes, welche die DickemeBeinrichtung erreicht, berechnet, d.h. in Abhingigkeit des Ende- 
erreicht-ldendif ikationssignals, welches durch die DickemeBeinrichtung erzeugt wird, genauer gesagt 
werden die vergangenen zeitsequentiellen Daten der Hitze, welche durch den Erhitzer erzeugt werden 
(in diesein Falle die Daten u(k-2), u(k-1) und u(k)), welche die durch die GroBe der in dem Speicher 

50 (104) gespeicherten Totzeit L bestimmt werden, der Zustandsvorhersageeinrichtung (106) zugef uhrt 

und der Zustandsveranderungsbetrag l(k+1) wird aus dem Eingang u(k) vom Zeitpunkt (WL) bis zum 
Zeitpunkt t* +1 erzeugt; 

(7) der Ausgang e AL [x,(k+1),©(k+1)] T des Zustandswechslers (105) und der Ausgang l(k+1) der Zu- 
standsvorhersageeinrichtung (106) werden in der Addiereinrichtung (107) addiert, welche den Zu- 

55 ^ A 

standsschfitzwert [X x (k+1), X(k+1)] T zum Zeitpunkt t^ erzeugt; auf diese Weise kann der Zustand- 

sch§tzwert zum Zeitpunkt t^ durch Integration in dem Zustandwechsler (105) und der 

Zustandsvorhersageeinrichtung (106) fur die Totzeit L erhalten werden, obwohl der Betriebsrechner 
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(103) nur den Zustandssch&tzwert zum Zeitpunkt infolge der Totzeit L erhalten kann; der EinfluB 
der PhasenverzSgerung infolge der Totzeit Lkann durch diese Operation eliminiert werden; 
(8) der Betrag u(k+1) der Hitze, welche durch den Erhitzer (12a) vom Zeitpunkt t^ bis zum nSchsten 
Berechnungszeitpunkt tk +2 erzeugt wird, wird durch die folgende Gleichung unter Verwendung derZu- 
standruckkopplungsverstarkung (ff, FJ bestimmt 



10 



u(k+1) = -fiX x (k+1)-F 2 X(k+1) (41) 

A A 

Die Addiereinrichtung(107)gibtdenZustandsschatzwert[x (k+1), x (k+1)] T zum Zeitpunkt tk +1 
an die Befehlseinrichtung (108) fur die durch den Erhitzer erzeugte Hitze ab; die Befehlseinrichtung 



15 
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(108) multipliziert den Zustandsschatzwert [ x (k+1),X(k+1)] T mit der Zustandsruckkopplungsverstar- 
kung zur Bestimmung eines Befehlswertes der durch den Erhitzer erzeugten Hitze; und 
(9) die obige Kontroilberechnung wird, nachdem der nSchste erfaSte Wert y(k+2) der Filmdicke von 
dem Abtaster (1 00) zum Zeitpunkt t=t^ 2 der Berechnungsausf uhrung erhalten wurde, ausgef uhrt, wenn 
die Dickeme&einrichtung entlang der Filmbreite nach Ablauf der Zeitperiode T bewegt wird und das ent- 
gegengesetze Ende erreicht, wobei die Gleichungen (29),(30),(33) t (34) I (37),(38) t (39) t gegeben sind, 
wie folgt: 

u (k*l) = ^ u> <k)*T l u(k-3)*T 2 u<k-2) ( 2 Q) 



u (kM)=u (k* 1 ) ♦KlydtM >-C<j (k*lH 



(30) 



25 



u (kM) = 0 u> <k)*r t u(k-2)*r 3 u(k-l) 



u> (kM) = u (kM)*K(y(kM)-Cu (k*l)3 



(33) 
(34) 
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KkU) = r ;t ^'j^"^nda u(k-2)^e AT I § 



at r 1a< 
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Mil <r u(k-l )♦ Bda u(k) 

o 



(37) 



KkM) = e ;T J^Bdau(k-l)»f pBdffji(k) (38) 
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X, (k*l) 



X<kH> 



X,(k»l) 



u (k*l ) 



♦ l(k*l) 



(39) 



so 8. Vonrichtung nach Anspruch 7, bei der die Ruckkopplungsmatrix F der Betriebsmengenbefehlseinrichtung 
(108) so ausgewahlt wird, daR alle Eigenwerte der Matrix (A-BF) in einem stabilen Bereich Negen. 

9. Vorrichtung nach Anspruch 4, bei der der Betriebsrechner (108) einen Sch§tzwert ©(k+1) fur©^) aus 
den folgenden Gleichungen berechnet: 
55 Die Berechnungsgleichung fur den SchStzwert ©(k+1) fur das Ende ® ist gegeben durch 
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T-U*Lc 

♦ I e A,| Bd7 u(k-2> (55) 
u <k*i) = u Ck*l) «■ K B Cy(k*l)-Ctj (k*l) ) (56) 

wobei 

(t >B = e A <T- L B + Lc) 

K B = die Verstarkungsmatrix des Betriebsrechners (1 03), und die Rechnungsgleichung fur den Schatzwert 
©(k+1) fur das Ende C gegeben ist durch 
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o <k+l) = 0 C C (k) «■ J™ C c Al Bd3/ u(k-l)_ 

♦ ?' L %Mat> u(k-2) (57) 

£ Ck*l> = cj <k*l) ♦ KcCy(k+l)-Cw (k>l)] (58) 
25 wobei *♦ 16 

*c = e A rr-L c + L B ) 
Kc= die Verstarungsmatrix des Betriebsrechners ist. 

10. Vorrichtung nach Anspruch 1, bei der die Kontrolleinrichtung die folgenden Bestandteile aufweist: 
so (1) der erfa&te Wert y(k+1) der Filmdicke (ein aus y^k+1, y 2 (k+1), y 3 (k+1), y 4 (k+1) und y 5 (k+1) beste- 

hender Vektor) wind durch die Dickeme&einrichtung (10) und den Abtaster (100) bei der Berechnungs- 
ausfuhrungszeitt=t k+1 desZeitintervalles Terhalten; der Abtaster (100) schlie&tbei jeder Berechnungs- 
ausfuhrungszeit t=t k+1f d.h. der Abtaster (100) schlie&t jedesmal, wenn die Dickeme&einrichtung (10) 
das Ende ® Oder ©des Filmes, wie im Figur 8 gezeigt, erreicht; wenn die Dickeme&einrichtung (10) 
35 das Ende ® Oder ©des Filmes erreicht, erzeugt die Dickeme&einrichtung (1 0) das Ende-erreicht Iden- 

difikationssignai d, welches das durch die Me&einrichtung erreichte Ende anzeigt; 
(2) ein Wert y 3 (k+1) des erfa&ten Filmdickewertes y(k+1) wird der Subrahiereinrichtung (101) zuge- 
fuhrt, welches die Dickenabweichung e(k+1) = r 3 (k+1)-y 3 (k+1 ) zwischen dem erfa&ten Wert y 3 (k+1) und 
einem Dickesollwert r 3 (k+1) erzeugt; 
40 (3) dem Integrator wird die Dickenabweichung e(k+1) aus der Subtrahierschaltung (101) zugefuhrtund 

er erzeugt einen zeitintegrierten Wert der Dickenabweichung aus der folgenden Gleichung; 

X,(k+1) = X l (k) + 0.5(t k . 1 >t k ){e(k) + 6(k + 1)) (69) 
wobei e(k) die Dickenabweichung bei der letzten Dickenerfassungszeit (t=y und X,(k) ein Ausgang des 
Integrators (102) bei t=tk ist; 

45 der Integrator (102) umfa&t eine Funktion eines Kompensators externer Stdrungen und dientzum Aus- 

gleich externer Hitze, welche die Dicke y 3 verandert, durch von dem Erhitzer erzeugte Hitze, so da& 
die Dicke y 3 immer bei einem Soll-Wert gehalten wird; 

(4) wenn eine Dickeme&einrichtung ein Ende des Filmes erreicht, erzeugt die Dickeme&einrichtung 
das Ende-erreicht Idendifikationssignal d. & (k+1) wird aus den Gleichungen (55) und (56) und (57) 

so und (56) in Reaktion auf das Idendifikationssignal d berechnet; die vergangene Datensequenz u(k-2) 

und u(k) der durch den Erhitzer erzeugten Hitze, welche in einem Speicher (104) zusammen mit dem 
erfa&ten Filmdickewert y(k+1) gespeichert ist, wird dem Betriebsrechner zugefuhrt welcher einen 
Schatzwert von X(t^ r L) = A(k+1) der Zustandsvariablen zum Zeitpunkt t( k+1 -L) erzeugt, welcher vor 
dem Zeitpunkt t k+1 urn die Totzeit L liegt, welche durch das Ende-erreicht Idendifikationssignal d be- 

55 stimmt wird, das durch die Dikkeme&einrichtung erzeugt wird; 

(5) bei der Berechnung des ersten Terms auf der rechten Seite der Gleichung_(66) wird der Zustands- 
schStzwert [x,(k+1),©(k+1)] T zum Zeitpunkt (WL) mit einem Koeffizienten zum Wechsel des Zu- 
stands bei der Durchschnittstotzeit C multipliziert, welche durch die Gleichung (68) bestimmt wird, urn 
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den Zustandsschatzwert e*-" [Xrfk+I), ©(k+1 )]T zum Zeitpunkt t^ zu erhalten; der Ausgang X,(k+1) des 
Integrators (102) und der Ausgang co(k+1) des Betriebsrechners (103) werden dem Zustandswechsler 
(1 05) zugef uhrt, welcher sie mit dem Koeff izienten zum Wechsel des Zustands bei der Durchschnitts- 
totzeit Cmultipliziert, um den Zustandsschatzwert zum Zeitpunkt W zu erhalten; die GroBe der Totzeit 
C nimmt den durch-schnitt lichen Wert der Totzeiten fur belde Enden des Filmes an, wie durch die Glei- 
chung (68) beschrieben ist; der Zustandwechsel durch den Eingang u(k), welcher in der Zeitdomane 
lediglich fur die Durchschnittstotzeit Cangelegt wird, wird durch den zweiten Term l(k+1) auf derrechten 
Seite der Gleichung (66) ausgedruckt und eine Korrektur wird daher durch die Zustandsvorhersage- 
einrichtung (106) ausgefuhrt; 

(6) der zweite Term f(k+1 ) auf der rechten Seite der Gleichung (66) druckt einen Bet rag des Zustands- 
wechsels fur die Zeiteingabedatensequenz u(k-1 ) und u(k) aus, welche der Zeitdomane fur die Durch- 
schnittstotzeit C von Zeitpunkt (WL) bis zum Zeitpunkt ^ angelegt werden; T(k+1 ) wird aus der Glei- 
chung (67) unterVerwendung der Durchschnittstotzeit E berechnet; die vergangene Datensequenz der 
durch den Erhitzer erzeugten Hitze (in diesem Fade zwei Daten u(k-1) und u(k)), welche durch die Grd- 
Be der in dem Speicher (104) gespeicherten Totzeit C bestimmt werden, werden der Zustandsvorher- 
sageeinrichtung (106) zugef uhrt und der Zustandsendebetrag l(k+1) durch den Eingang u(k) vom Zeit- 
punkt (WR bis zum Zeitpunkt. 

(7) der Ausgang e*4X,(k+1) t ©(k+1)] T des Zustandswechslers (105) und der Ausgang l(k+1 ) Zustands- 
vorhersageeinrichtung werden in der Addiereinrichtung (107) addiert, welche den Zustandsschatzwert 

^\ 

[X i (k+1), X(k+1)] T zum Zeitpunkt t^ erzeut; obwohl der Betriebsrechner (103) lediglich den Zustands- 
schatzwert zum Zeitpunkt WL infolge der Totzeit t erhalten kann, kann der Zustandsschatzwert zum 
Zeitpunkt t k+1 durch Integration in dem Zustandswechsler (105) und der Zustandsvorhersageeinrich- 
tung (106) fur die Totzeit L erhalten werden; der EinfluB der Phasenverzogerung infolge der Totzeit L 
kann durch diese Operation eliminiert werden; 

(8) ein Betrag u(k+1) der durch den Erhitzer erzeugten Hitze vom Zeitpunkt t^ bis hin zum nachsten 
Berechnungszeitpunkt t*+ 2 wird durch die fblgende Gleichung unter Verwendung der Zustandsruck- 
kopplungsverstirkung (f 1f F 2 ) bestimmt: 

u(k+1) = -f,§ x (k+1)-F 2 X(k+1) (41) 

die Addiereinrichtung (107) gibt den Zustandsschatzwert [ x (k+1), X(k+1)] T zum Zeitpunkt t^ an die 
Befehlseinrichtung (108) fur die durch den Erhitzer erzeugte Hitze ab; die Befehlseinrichtung (108) 

A A 

multipliziert den Zustandsschatzwert [ x (k+1), X(k+1)] T mit der Zustandsruckkopplungsverstarkung 
zur Bestimmung eines Befehlswertes der durch den Erhitzer erzeugten Hitze; und 

(9) die obige Berechnung wird ausgefuhrt nachdem der nachste erfaBte Wert y(k+2) der Filmdicke von 
dem Abtaster (100) zum Zeitpunkt t=tk* 2 der Berechnungsausfuhrung erhalten wird, wenn die 
DickemeBeinrichtung nach der Zeitperiode T entlang der Breite des Filmes bewegt wird und das ent- 
gegengesetzte Filmende erreicht, wobei die Gleichungen (55), (56),(57), (58), (66), (67),(68), gegeben 
sind folgt: 



<k*l) = 0 B C (k) * I e e*"Bd?u(k-l) , S T ' L e X Hdv u(K-2) (55) 

fl)|] 

u(k^l) = ui • k ♦ ) ) ♦ K a ly<k*l )-Co7 (k*L ) ) (55) 
wobei 

4> B = eA(T-LB + Lc) 
Ke= die Verstfirkungsnatrix des Betriebsrechners ist. 
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Hie 



u(k*l> = 0 rju (k) ♦ J g*\\Mv u<k-l ) 



r T-LciI.n 
♦ J c A ' I3d ?f u(k-2 > 



u(kM) = u> <k*l) ♦ K c [y(k>l >-Cw (k*l)] 



(57) 
(58) 



fo wobei 

*c = e*(T-L c + L B ) 
Kc= die Verstarkungsmatrix des Betriebsrechners ist. 
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X, (k*L) 



X(k+1 > 



X, (k*l> 



us (k*l) 



♦ l(k*l) 



(66) 



as T(k*l) = e**J* e A ' f 6do u(k-l) 



J T e Rr Gd o 

O 



ii (k) 



(67) 
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L = (Lb + L c )/2 (68) 



35 Revendications 

1. Dispositif de commande de I'epaisseur d'une couche qui utilise des Equations d'6tat pour commander un 
appareii de mouiage par extrusion ou un appareil de moulage par 6coulement de couches comportant 
une f iliere ayant un mecanisme qui commande, dans la largeur de la couche, un debit de plastique fondu, 

40 et une jauge d'epaisseur pour detector une variation d'epaisseur de la couche apres I'ecoulement d'un 

temps mort Li correspondent au temps necessaire au mouvement de la couche entre la f iliere et la jauge 
d'epaisseur, comprenant un soustracteur (101) pour produire une difference entre une vaieur d'epaisseur 
detectee par la jauge d'epaisseur (10) dans une position predeterminee dans la largeur de la couche et 
une vaieur de consigne d'epaisseur dans la position predeterminee, ou ladite jauge d'epaisseur se de- 

45 place dans un mouvement alternatif dans la largeur de la couche, caracterise par un integrateur (102) 

pour integrer dans ie temps la difference d'epaisseur produite par ledit soustracteur (1), une memoire 
(104) pour stocker des donnees de sequence de temps passe de quantites de fonctionnement d'un dis- 
positif de reglage d'epaisseur (12b) pendant un temps mort L egal a la somme du temps mort et d'un 
temps L 2 jusqu'a ce que la jauge d'epaisseur (10) atteigne une extremite de la couche apres detection 

so de I'epaisseur a la position predeterminee, un calculates operationnel (103) pour produire les donnees 

de sequence de temps passe de quantites de fonctionnement du dispositif de reglage d'epaisseur (12b) 
stockees dans ladite memoire, et une vaieur estimee de variable d'etat a un temps precedant celui ou la 
vaieur de consigne de la vaieur detectee de I'epaisseur de la couche a ete introduite du temps mort L, un 
modificateur d'etat (105) pour recevoir une sortie dudit integrateur (102) et une sortie dudit calculates 

55 operationnel (103) et multiplier chacune desdites sorties dudit integrateur (102) et calcu I ateur operationnel 

(103) par un coefficient pour modifier I'etat du temps mort L pour produire une vaieur estimee d'etat a un 
temps predetermine, un dispositif de prevision d'etat (106) pour recevoir les donnees de sequence de 
temps passe de quantites de fonctionnement du dispositif de reglage d'epaisseur (12b) stockees dans 



37 



EP 0 329 157 B1 



ladite memoire (104) pour produire une variation d'etat basee sur I'etablissement de Pentree, depuis un 
certain temps jusqu'a un temps apres I'ecoulement du temps mort L, un additionneur (107) pour addition- 
ner une sortie dudit modif icateur d'etat (105) et une sortie dudit dispositif de prevision d'etat (106) pour 
produire la valeur estimee d'etat au temps predetermine, et un dispositif de commande de quantite de 
fonctionnement (108) pour multiplier une valeur estimee d'etat a un certain temps produite par (edit ad- 
ditionneur par un gain de retroaction d'etat pour produire une valeur de commande de quantite de fonc- 
tionnement dudit dispositif de reglage d'epaisseur (12b). 

Dispositif de commande selon ia revendication 1 , ou lesdites equations d'etat sont donnees par 

x(t)=Ax(t) + Bu(t) (2) 
y(t) = Cx(t-L) (4) 

ou x est un vecteur d'etat, u est un vecteur d'entree dans lequel u(t) = [u^t), u 2 (t), u 3 (t), u 4 (t), u 5 (t)] T (ou T 
represente transposition), y est un vecteur de sortie dans lequel y(t)= [y^t), y 2 (t), y 3 (t), y 4 (t), y 6 (t)] T , et les 
equations d'etat (2) et (3) sont contrfilables et observables. 

Dispositif de commande selon la revendication 2, ou ledit integrates (102) realise le calcul suivant 



X,(t) = - J*CiX(T )dr - j'?iX(r )dT (6) 
ou Ci est la i-ieme rangee de la matrice C. 

Dispositif de commande selon la Revendication 3, ou lorsque c>(t)=x(t-L), ledit calculateur operationnel 
(103) calcule une valeur estimee X(t-L) pour x(t - L) sur la base des equations pour ©(t) suivantes: 

o>(t)= A©(t) + Bu(t-L) (18) 
y(t) = C<D(t) (19) 

Dispositif de commande selon la revendication 4, oil ladite equation (18) est I'equation discrete suivante: 
o> (t K -i) = e A (t k .!-t k )w (t K ) ♦ J*e A ( t K -i-T > 



Bu( T -L)dr 

Dispositif de commande selon la revendication 6, ou lorsque r| = k+i-T, une valeur estimee ©(k+1) pour 
©(t k+1 ) est calculee par les equations suivantes: 
dans le cas de 2^L<T, 

cj (k*l)=0 Z (k) + T x u(k-3) + r a u(k-2) (29) 
w (k+l) = £> (k+1 )+Kty(k*l)-Cw (k>l)] (30) 

ou K est une matrice de gain de retroaction dudit calculateur operation nel (108), et dans le cas de T<L<2T, 
Z (k+l) = 0 (2 (k) + T x u(k-2) + r a u(k-l) (33) 

£ Ck+l) = w (k+1 )+KCy(k+l>-Cu (k+1)] (34) 

ou 

((> = eAT (23) 

r i = x% A ^Bd5? 

r 2 = I™e A7 »Bd7 



(24) 
(25) 
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m = 3T - L (26) 

Dispositif de commande selon revendication 1 , ou ledit dispositif de commande comporte les constituants 
suivants : 

(1) une valeur d6tect6e y(k+1) d'epaisseur de couche (vecteur constitu6 de yi(k+1), y 2 (k+1), y 3 (k+1), 
y 4 (k+1) et y 5 (k+1)) est obtenue au travers de la jauge d'epaisseur (10) et un 6chantillonneur (100) au 
temps d'ex6cution de calcul t=t k+1 de I'intervalle de temps T; l'6chantillonneur (100) se ferme pour cha- 
que temps d'ex6cution de calcul t=t k+1 , c'est-a-dire, rechantillonneur (100) se ferme chaque fois que 
la jauge d'epaisseur (1 0) atteint I'extremite ® ou © de la couche representee dans la Figure 8 ; de 
plus, lorsque la jauge d'epaisseur (10) atteint I'extremite ® ou © de la couche, la jauge (10) produit 
le signal d'identification d'extremite d'arriv6e d qui indique l'extr6mite que la jauge d'epaisseur (10) a 
atteinte; 

(2) une valeur y 3 (k+1 ) de la valeur y(k+1) d6tect6e d'epaisseur de couche est fournie au soustracteur 
(101) qui produit la d6viation d'6paisseur e(k+1)=r 3 (k+1)-y 3 (k+1) entre la valeur detectee y 3 (k+1) et une 
valeur de consigne d'epaisseur r 3 (k+1); 

(3) I'int6grateur (102) recoit la deviation d'epaisseur e(k+1) du soustracteur (101) et produit une valeur 
int6gr6e dans le temps de la deviation d'epaisseur, £ partir de liquation suivante: 

X,(k + 1) = X,(k) + 0.5& ♦ ^ - t0{e(k) ♦ e(k + 1)} (40) 
ou e(k) est ia deviation d'epaisseur au dernier temps de detection d'epaisseur (t=t k ) et X^k) est une 
sortie de I'integrateur (102) & t=t k ; 

I'int6grateur (102) inclut une fonction d'un compensateur de perturbation externe, et sert £ 
compenser de la chaleur externe faisant varier l'6paisseur y 3 par de ia chaleur g6ner£e par ('element 
de chauffage, si bien que repaisseur y 3 est toujours ajustee pour etre une valeur de consigne; 

(4) lorsque la jauge d'epaisseur atteint I 'une ou I'autre des extr6mites de la couche, la jauge d'epaisseur 
produit le signal d'identification d'extremite d'arrivee d; ©(k+1) est calcuie £ partir des equations (29) 
et (30) ou (33) et (34) en r6ponse au signal d'identification d; plus particulierement, dans les equations 
(29) et (30) pour les donn6es de la sequence de temps passe de chaleur generee par i'6l6ment de 
chauffage (108), stockees en memoire (104), u(k-3) et u(k-2) sontfournies au calculateur op6rationnel 
(1 03), tandis que dans les equations (33) et (34), u(k-2) et u(k- 1 ) ainsi que ia valeur detectee d'epaisseur 
de couche y(k+1 ) sont fournies au calculateur op6rationnel (1 03), qui produit une valeur estim6e X(ttc +1 - 
L)=o)(k+1 ) de la variable d'etat au temps t( k+1 - L) pr6c6dant le temps t k+1 du temps mort L determin6e 
par le signal d'identification d'extremite d'arrivee d produit par la jauge d'epaisseur; 

(5) dans le calcul du premier terme du membre droit de ('equation (39), la valeur est i m6e d'etat [Xj((k+1 ), 
©(k+1)] T au temps t( k+1 -L) est multiplies par un coefficient e^ pour modifier l'6tat au temps mort L, pour 
obtenir la valeur estim6e d'etat e^X^k+l^k+l)] 1 " au temps tk +1 ; c'est-£-dire la sortie X,(k+1) de I'in- 
tegrateur (1 02) et la sortie ©(k+1 ) du calculateur op6rationnel (1 03) sont fournies au modif icateur d'etat 
(105) qui les multiplie par le coefficient, pour changer I'etat au temps mort L, determine par le signal 
d'identification d'extremite d'arrivee d de la jauge d'epaisseur, pour obtenir la valeur estim6e d'etat au 
temps t k+1 ; etant donn6 que la grandeur mort L diff6re selon I'extr6mit6 de la couche que la jauge 
d'epaisseur atteint, le coefficient e^ diff6re selon la position de la jauge d'epaisseur £ l'ex6cution du 
calcul, c'est-a-dire, le signal d'identification d'extremite d'arrivee d de la jauge d'epaisseur; 

le changement d'6tat par I'entr6e u(k) applique en domaine temps uniquement pendant le temps 
mort L est exprim6 par le second terme l(k+1 ) du membre droit de l'6quation (39) et la correction corres- 
pondante est r6aiisee par le dispositif de prevision d'etat (106); 

(6) le second terme l(k+1) du membre droit de l'6quation (39) exprime une quantite de modification 
d'6tat pour les donn6es d'entr6e de sequence de temps u(k-2), u(k-1), et u(k) appliqu6es au domaine 
de temps, du temps (t^ - L) au temps t k+1 : l(k+1) est calcuie £ partir de l'6quation (37) ou (38), en fonc- 
tion de I'extremite de la couche que la jauge d6paisseur atteint, c'est-d-dire, en fonction du signal 
d'identification d'extremite d'arrivee produit par la jauge d'epaisseur, plus particulierement, les don- 
nees pass6es, s6quentielles selon le temps, de la chaleur generee par reiement de chauffage (dans 
ce cas, trois donn6es de u(k-2), u(k-1), et u(k)), d6termin6es par la grandeur du temps mort L stock6 
dans la memoire (104), sontfournies au dispositif de prevision d'etat (106), et la quantite de variation 
d'etat l(k+1) par I'entr6e u(k) du temps (t^ - L) au temps t k+1 est produite; 

(7) la sortie eAL[X,(k+1 ),w(k+1 ) T du modif icateur d'etat (1 05) et ia sortie l(k+1 ) du dispositif de prevision 
d'etat (106) sont additionn6es dans I'additionneur (107) qui produit la valeur estim6e d'etat 

IX i (k+1),X(k+1)] T au temps tfc +1 ; de cette maniere, bien que le calculateur operationnel (103) puisse 
n'obtenir que la valeur estim6e d'etat au temps tk +1 -L £ cause du temps mort L, la valeur estim6e d'etat 
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au temps t k+1 peut dtre obtenue par integration dans le modif icateur d'etat (105) et le dispositif de pre- 
vision d'etat (106) pour le temps mort L; ('influence du retard de phase du au temps mort L peut etre 
eliminee par cette operation; 

(8) une quantite u(k-1) de chaleur generee par I'element de chauffage (12a) du temps t k+1 au prochaln 
temps de calcul t k+2 est definie par I'equation suivante qui utilise le gain de retroaction d'etat (f 1t FJ; 

u(k+1) = -f 1 f I (k+1)-F 2 X(k+1) (41) 

A a 

I'additionneur (107) fournit la valeur estimee d'etat [ x (k+1),X(k+1)] T au temps t^a un dispositif de 
commande (108) pour de la chaleur generee par I'element de chauffage; le dispositif de commande 

a A 

(108) multiplie la valeur estimee d'etat [ x (k+1),X(k+1)] T par le gain de retroaction d'etat, pour deter- 
miner une valeur de commande de chaleur generee par I'element de chauffage; et 

(9) le calcul de commande ci-dessus est execute apres que la valeur detectee suivante y(k+2) d'epais- 
seur de couche soit obtenue de I'echantillonneur (100) au temps t=t k+2 de I'execution du calcul, lorsque 
la jauge d'epaisseur est deplacee dans la largeur de la couche apres la periode de temps T, et qu'elle 
atteint I'extremite oppos6e de la couche, lesdites equations (37), (38), (39) etant les suivantes : 



KK M) = c*« aT '^-?*'ijda u(k-2)*e* T J p r 

lW ii ( k- I ) ♦ £e**Bdff u(k> (3 7) 



1 (k+l) 


- e Sr C 


~e A<r Bda u(k-L)*J |**Bda u(k) 


(38) 


X t (k*l ) 




X,(k*l) N 






X(k*l) 




w (k + 1) 


♦ I(k*l) 


(39) 



Dispositif de commande selon la revendication 7, dans lequel la matrice de retroaction F dudit dispositif 
de commande de valeur de fonctionnement (1 08), est choisie de facon & ce que chacune des valeurs pro- 
pres de la matrices (A-BF) soient dans une region stable. 

Un dispositif de commande selon la revendication 4, dans leque) (edit calculateur operationnel (108) cal- 
cule la valeur estimee <&(k+1 ) pour o>(t k+1 ) & partir des equations suivantes: 

I'equation de calcul de la valeur estimee ©(k+1 ) pour I'extremite ® est donnee par 

Z (k*l) = * »u (k) ♦ £e A, 'Bd* u(k-l) 
T-L 8 *Lc 

♦ J e An Bd* u(k-2) (55) 
Rig 



u (k+1) = u ♦ K B ty(k+l)-CuT (k+1)] (56) 

ou 
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<fe = eA<T-LB + Lc) 
K B = matrice de gain du calculateur operationnel (103), 

liquation de calcul de la valeur estimee <S(k+1 ) pour I'ext remite © est donnee par 
U (k+l) = 0 c u (k) ♦ J 0 e A "Bd7 u(k-l) 

♦ S T ~ L t*k%dy u(k-2) (57) 

me 



u> (k+l) = cj (k+l) ♦ KcCyCk+l)-Cu (k+l)3 (58) 



+0 = 6^-10 + 18) 

Ko= matrice de gain du calculateur operationnel. 

Dispositif de commande selon la revendication 1, ou ledit dispositif de commande comprend les consi- 
tuants suivants: 

(1) la valeur detectee d'epaisseur de couche y(k+1) (vecteur constitue de y n (k+1, y 2 (k+1), y 3 (k+1), 
y 4 (k+1)et y 5 (k+1)), est obtenu au travers de la jauge d'epaisseur (10) et I'echantillonneur (100), au 
temps d'execution du calcul t=t k+1 de I'intervalle de temps T; rechantillonneur (1 00) se ferme pour cha- 
que temps d'execution du calcul t=t k+1 , c'est-a-dire, rechantillonneur (100) se ferme chaque fois que 
la jauge d'epaisseur (1 0) atteint I'extremite ® ou ©de la couche representee dans la Figure 8 ; lorsque 
la jauge d'epaisseur (10) atteint I'extremit6 ® ou ©de ia couche, la jauge (10) produit le signal d' iden- 
tification d'extremite d'arriv6e d qui indique I'extremite atteinte par la jauge; 

(2) une valeur y 3 (k+1 ) de la valeur detectee d'epaisseur de couche y(k+1 ) est fournie au soustracteur 
(101) qui produit la deviation d'epaisseur e(k+1) = r 3 (k+1)-y 3 (k+1) entre la valeur d6tect6e y 3 (k+1) et 
une valeur d'epaisseur de consigne r 3 (k+1); 

(3) I'integrateur (1 02) recoit la deviation d'epaisseur e(k+1 ) du soustracteur (101) et produit une valeur 
integree dans ie temps de la deviation d'epaisseur a partir de requation suivante: 

X l (k+1) = X l (k) + 0.5(t k . 1 -t k ){e(k) + 8(k+1)} (69) 
ou e(k) est la deviation d'epaisseur au dernier temps de detection d'epaisseur (t=t0 et X|(k) est une 
sortie de I'integrateur (102) a t=t k ; 

I'integrateur (102) inclut une fo notion d'un compensateur de perturbation externe et sert a 
compenser de la chaleur externe faisant varier l'6paisseur y 3 par de la chaleur g6n6r6e par t'6l6ment 
de chauffage de facon a ce que I'epaisseur y 3 soit toujours ajustee pour etre une valeur de consigne; 

(4) lorsque la jauge d'epaisseur atteint Tune ou I'autre des extremites de la couche, la jauge d'epaisseur 
produit le signal d'identification d'extremite d'arrivee d; o(k+1) est calcuie a partir des equations (55) 
et (56) ou (57) et (58) en reponse au signal d'identification d; les donn6es de sequence de temps passe 
u(k-2) et u(k-1 ) de chaleur gen6r6e par l'6l 6ment de chauffage, stock6es en m6moire (1 04) avec la va- 
leur detectee d'epaisseur de couche y(k+1), sontfournies au calculateur operationnel, qui produit une 
valeur estimee X(t k+1 -L)=m(k+1 ) de la variable d'etat au temps t( k+1 -L) pr6c6dant le temps t^ du temps 
mort L determine par le signal d'identification d'extremite d'arrivee d produit par la jauge d'epaisseur; 

(5) dans le calcul du premier terme du membre droit de requation (66), la valeur estimee d'etat 
[X,(k+1 ),ca(k+1 )] T au temps t( k+1 -L) est multiple par un coefficient e^ pour modifier I'etat par le temps 
mort moyen C d6f ini par requation (68) pour obtenir la valeur estimee d'etat e AL [X,(k+1),©(k+1)] T au 
temps t k+1 ; la sortie X,(k+1) de I'integrateur (102) et la sortie u(k+1) du calculateur operationnel (103) 
sont fournies au modificateur d'etat (105) qui les multiplie par le coefficient pour modifier i'etat par le 
temps mort moyen L pour obtenir la valeur estim6e d'6tat au temps t^, la grandeur au temps mort C 
adopte la valeur moyenne des temps morts pour les deux extremites de la couche, comma d6crit par 
requation (68); 

la modification d'etat par I'entr6e u(k) appliqu6e en domaine de temps uniquement pendant le 
temps mort moyen L, est exprim6e par le second terme T(k+1) du membre droit de requation (66) et la 
correction correspondante estfaite par le dispositif de prevision d'etat (106); 

(6) le second terme T(k+1) du membre droit de l'6quation (66) exprime une quantite de modification 
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d'etats pour les don nees d'entree de sequence de temps u(k- 1 ) et u(k) appliquees au domaine de temps 
du temps mort moyen C, du temps (WQ au temps t k+1 : T(k+1) est calcule a partir de I'equation (67) 
en utilisant le temps mort moyen C; les donnees de sequence de temps passe de la chaleur generee 
par I'element de chauffage (dans ce cas, deux donnees de u(k-1) et u(k)) determinees par la grandeur 
du temps mort Cstocke dans la memoire (104), sontfournies au dispositif de prevision d'etat (106) et 
la quantite devariation d'etat l(k+1) par I'entree u(k) du temps (WC) au temps t^; 

(7) la sortie e^[X,(k+1),©(k+1)] T du modificateur d'etat (105) etla sortie l(k+1) du dispositif de prevision 

d'etat, sont additionnees dans I'additionneur (107) qui produit la valeur estimee d'etat [Xi(k+1), 
Xfk+I)] 1 au temps T k+1 ; bien que le calculates operationnel (103) ne puisse obtenir que la valeur es- 
timee d'etat au temps T k+1 -L en raison du temps mort L; la valeur estimee d'etat au temps T k+1 peut 
Stre obtenue par integration dans le modif icateur d'etat (105) et le dispositif de prevision d'etat (106) 
pour le temps mort C; I'influence du retard de phase dO au temps mort L peut 6tre eliminee par cette 
operation; 

(8) une quantite u(k+1) de chaleur generee par I'element de chauffage du temps t k+1 au temps de calcul 
suivant t k+2 est definie par I'equation suivante, qui utilise le gain de retroaction d'etat (f lf FJ: 

u(k+1) = -f^ x (k+1)-F 2 X(k + 1) (41) 

a A 

I'additionneur (1 07) fournit la valeur estimee d'etat [ x (k+1),X(k+1)] T au temps t k+1 a un dispositif de 
commande (108) pour de la chaleur generee par I'element de chauffage; le dispositif de commande 

a A 

(108) multiplie la valeur estimee d'etat [ x (k+1),X(k+1)] T par le gain de retroaction d'etat pour deter- 
miner une valeur de commande de chaleur generee par I'element de chauffage; et 

(9) le calcul ckdessus est execute apres que la valeur detectee suivante y(k+2) d'epaisseur de couche 
est obtenue de I'echantillonneur (100) au temps t=tfc +2 d'execution de calcul, lorsque la jauge d'epais- 
seur est dSplacee dans la largeur de la couche apres la periode de temps T et atteint I'extremite op- 
posee de la couche, lesdits equations (55), (56), (57), (58), (66), (67), (68) etant les suivantes : 



30 



u Ck*l> = 0 B u (k) ♦ J™e* ? »Bdv u(k-l) ♦ j T L °^Bd v u(k-2) (55) 



u(k+l> = c5<k*l) ♦ K B Cy<k*l>-Co> <k*l>] 



(50) 



ou 
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ou 



* B = eA(T-L B + Lc) 

u(kM) = 0 c w(k) ♦ £ c c A, Mv u*k-'. ) 
♦ J^e^Bd* u(k-2) 
£ <k+l) = w(kM) ♦ K c [y<k^!)-Co7(k^l)] 

<|>c = eA(T-L c + L B ) 



(f)7) 



<58) 





Xi(k*l> 
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X, 
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(k*l) 
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♦ I(k+i> 



(66) 
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♦ J% x/ Bdcr u(k) 

C = (Lb + Lc)/2 (68) 
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